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SUMMARY 


l 

’  A  multipulse  laser-schlieren  technique  was  used  to  study  the  time 
behaviour  of  instabilities  in  the  translational  front  of  planar  shock  waves 
in  ionizing  argon  in  the  10  x  18  cm  hypervelocity  shock  tube  at  UTIAS.  Using 
a  ruby  laser  in  a  combined  normal-lasing/giant-pulse  mode,  at  least  two 
monochromatic  light  pulses  were  obtained  about  7  apart.  As  well,  a 
schlieren  system  compatible  with  the  23  cm  dia  Mach-Zehnder  interferometer 
was  created  with  maximum  sensitivity  for  some  low-density  conditions.  The 
almost-sinusoidal  perturbations  were  found  to  change  with  time;  some  like  a 
standing  wave,  and  others  like  a  travelling  wave.  It  was  not,  however,  possible 
to  associate  the  observed  perturbations  with  either  form  of  wave  motion  in 
nearly  half  of  the  cases.  A  simple  graph  provides  a  way  to  predict  the  shape 
of  the  perturbations  from  the  shock-wave  parameters. 
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1. 


INTRODUCTION 


1.1  Significance  of  Planar  Flow  in  a  Shock  Tube 

A  simplifying  assumption  in  shock-tube  flows  is  that  the  flov  is  one¬ 
dimensional.  Except  for  thin  boundary  layers  in  tubes  of  large  cross-section, 
this  supposition  is  quite  reasonable  and  has  been  verified  experimentally  by 
many  researchers.  The  shock  front  under  such  circumstances  is  found  to  be 
planar  at  a  few  diameters  from  the  diaphragm  station.  Notable  exceptions  to 
this  expected  planarity  are  the  ionizing  shock  waves  in  argon  and  krypton 
produced  in  the  UTIAS  hypervelocity  shock  tube  (Ref.  1,2,3).  Interferograms  of 
such  shocks  display  a  translational  front  with  decidedly  sinusoidal  oscillations 
along  with  coupled  disturbances  in  the  subsequent  relaxation  regions.  The 
existence  of  this  phenomenon  seems  to  contradict  the  expected  shock-wave 
stability  in  a  shock  tube.  It  also  hinders  the  reduction  of  shock-structure 
data  from  the  interferograms.  It  has  been  found  that  the  addition  of  a  small 
amount  of  hydrogen  to  the  test  gas  renders  the  shock  front  and  the  flow 
planar  (Ref.  1,2,3).  To  date  no  satisfactory  explanation  has  been  put  forward 
for  this  anomaly  but  several  suggestions  have  been  made  by  previous  researchers 
pertaining  to  the  source  of  the  non-planarities,  mechanisms  of  propagation,  and 
recommended  investigations .  The  remaining  part  of  this  introduction  to  the 
present  study  will  l)  provide  a  precise  description  of  the  problem  and  relate 
the  significance  of  the  investigated  phenomenon;  2)  summarize  the  findings  and 
suggestions  of  previous  researchers  and  3)  state  the  purpose  of  the  present 
study  as  influenced  by  these  proposals. 

Shock  tubes  have  been  used  for  many  years  in  conjunction  with  various 
optical  devices  in  investigations  of  high-speed  shock  waves  and  high-temperature 
gas  flows  (Ref.  U).  The  monatomic  gas  argon  has  been  used  extensively  by 
researchers  studying  plasma  behaviour  and  ionization  effects  because  it  is 
readily  available  in  high  purity ,  safe  to  work  with,  free  of  dissociation  effects, 
relatively  easy  to  ionize,  and  has  thermodynamic  properties  that  are  well  known 
over  a  wide  range  of  temperatures  and  pressures.  Because  of  its  high  atomic 
weight  (Uo  amu)  its  sound  speed  is  not  high  320  m/s  at  room  temperature),  so 
that  higher  flow  Mach  numbers  can  be  obtained  for  the  same  driver  and  channel 
pressures  (Ref.  U).  Lastly,  in  an  argon  plasma,  the  vast  difference  between 
me  and  m^.  (m^.  a*  70,000  me)  constitutes  a  significant  kinetic  barrier  between 
the  electrons  and  the  argon  ions  and  atoms ,  making  possible  two  distinct 
temperatures  Te  and  T^. ,  even  though  the  species  are  mixed  (Ref.  1-3,  5-7). 
Researchers  utilizing  the  UTIAS  10  x  18  cm  combustion-driven  shock  tube  and  the 
23-cm  diameter  Mach-Zehnder  interferometer  have  made  use  of  some  of  these 
features  of  argon  to  study  supersonic  corner  expansion  of  plasma  (Ref.  8), 
refractive  index  of  ions  (Ref.  l),  shock  structure  (Ref.  1-3,  9-ll)»  shock- 
induced  boundary  layers  (Ref.  11-15)  and  oblique  shock-wave  reflections 
(Ref.  16-18). 

Except  for  the  thin  sidewall  boundary  layer,  it  is  usually  assumed  in  a 
theoretical  analysis  that  the  flow  in  the  unobstructed  channel  is  one¬ 
dimensional,  that  is,  all  state  variables  of  the  test  gas  at  every  point 
(x,y,z)  in  the  channel  are  functions  of  the  space  variable  x,  only  (Fig.  1). 

This  assumption  greatly  simplifies  the  mathematical  modeling  of  the  gas- 
dynamics  of  the  shock,  which  otherwise  would  be  too  unwieldy  or  would 
involve  prohibitive  computational  costs. 
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The  laboratory  coordinate  system  generally  used  is  shown  in  Fig.  1,  where 
the  (x,y,z)-axes  are  at  the  diaphragm  station  and  oriented  in  the  following  way 
The  x-axis  is  parallel  to  the  longitudinal  axis  of  the  shock  tube;  the  y-axis 
is  in  the  vertical  direction;  the  z-axis  is  the  third  axis  of  the  right- 
handed  coordinate  system.  The  shock-fixed  coordinates,  (x' ,y',z')  are  related 
to  the  lab  coordinates  by 

x'  =  x  -  Ut,  y '  =  y,  and  z '  =  z, 

where  U  is  the  shock  speed  and  the  shock  front  is  stationary  at  x'  =  x  -  Ut  =  0 
and  the  gas  travels  through  the  shock  in  the  negative  x'-direction. 

The  interferometer  receives  from  the  light  source  light  travelling  in 
the  z-direction,  and  its  field  of  view  represents  the  sum  of  all  (x,y)-planes 
over  the  shock-tube  width  of  10. l6  cm.  If  the  flow  is  one-dimensional,  each 
state  variable  is  constant  on  each  (y,z)-plane;  hence,  the  equivalent  term 
"planar"  flow.  Each  (y,z)-plane  appears  at  the  image  of  the  test  section  as 
a  straight  vertical  line,  so  the  presence  of  planar  flow  in  the  test  section 
implies  that  every  observable  variable  is  constant  along  such  vertical  lines 
in  the  resultant  picture.  Great  care  is  taken  in  the  alignment  of  the  inter¬ 
ferometer  to  ensure  that  photographic  light  from  the  light  source  passing 
through  the  test  section  is  in  fact  parallel  to  the  z-direction,  that  is, 
perpendicular  to  the  window  surfaces.  This  alignment  is  routinely  checked 
by  observing  that  light  from  the  interferometer's  point-light  source  returns  to 
the  original  point  after  partial  reflection  off  the  various  test  section  window 
surfaces.  A  discussion  of  errors,  tolerances  and  adjustments  of  the  optics  of 
the  interferometer  is  found  in  Refs.  19  and  20. 

It  follows  from  the  above  discussion  that  the  surface  of  a  shock  wave  in 
on-dimensional  flow  will  necessarily  be  planar,  causing  the  shock  to  appear  as 
a  very  thin  straight  vertical  line,  if  the  light  pulse  used  to  obtain  the 
picture  is  short  enough  to  "freeze"  the  motion.  This  line  is  readily  seen  in 
schlieren  shadowgraph  and  inter ferogram  (Fig.  2)  pictures,  but  must  also  be 
visualized  on  an  interferogram  by  connecting  a  series  of  discontinuous  fringe 
shofts  (Figs.  3-6).  The  abruptness  of  the  fringe  shifts  attests  to  the  very 
thin  shock  surface  and  shows  that  schlieren  and  shadowgraph  images  of  the 
translational  shock  do  not  reflect  the  actual  thickness  which  is  on  the  order 
of  10-5  cm,  or  a  few  mean  free  paths. 

There  are  two  ways  in  which  the  shock  surface  can  be  distinguished  as 
non-planar.  The  first  occurs  when  the  shock  surface  position  depends  upon  x 
and  y  but  not  on  z;  the  image  consists  of  a  curved  thin  line  and  represents  a 
reflected  translational  shock  front  (Ref.  l8)(C  in  Fig.  2).  The  second  non¬ 
planarity  occurs  when  the  shock  surface  is  three-dimensional,  when  the  position 
depends  on  z  as  well  as  x  and  y.  This  effect  manifests  itself  in  stray  shocks 
branching  out  from  the  main  image  or  in  a  distinctive  "twisted  ribbon" 
appearance  where  it  seems  that  the  observer  is  looking  edge-on  at  a  flexible 
twisted  sheet  (Fig.  3).  Recognition  of  three-dimensional  phenomena  from  inter- 
ferograms  (schlieren  or  shadow  photographs)  is  relatively  easy,  as  small 
deviations  from  a  thin  shock  front  can  be  observed. 

Since  the  flow  behind  the  shock  front  has  gone  through  the  wave,  it  can 
be  expected  that  the  flov  is  coupled  with  it.  The  flow  behind  a  planar  shock 
front  is  usually  planar.  This  is  easily  depicted  by  interf erograms ,  in  which 
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fringe  shifts  are  constant  along  vertical  lines  everywhere  in  the  field  of 
view.  Planar  flow  is  most  vividly  illustrated  when  the  no-flow  fringes  are 
rotated  to  he  horizontal,  that  is,  parallel  to  the  x'  direction.  In  this  case, 
every  fringe  is  congruent,  even  in  chemically- reacting  flows  at  high  temperature 
(Fig.  4).  In  plotting  a  profile  of  some  state  variable  through  a  shock-wave 
transition  (e.g.,  density  p  vs  x')  from  a  picture  such  as  Fig.  4,  one  can  be 
confident  that  at  every  position  x'  the  variable  had  that  value  over  the  whole 
(y' ,  z*) -plane.  That  is,  the  profile  applies  to  the  whole  volume  of  the  test 
section.  Similarly,  when  a  variable  such  as  relaxation  length  xg  is  quoted, 
there  is  no  ambiguity  because  every  fringe  provides  the  same  answer.  In 
addition,  one  may  justifiably  apply  to  the  observed  physical  phenomena  one  of 
the  previously  mentioned  theoretical  analyses  which,  even  though  one-dimensional, 
in  many  cases  is  very  complex  and  computationally  expensive.  Hence  interfero¬ 
metric  verification  of  planar  flow  in  the  shock  tube  is  essential  to  the 
reliable  interpretation  of  measured  data. 

1.2  Review 


While  studying  ionizing  argon  gas  flows  in  this  facility,  Bristow  (Ref.  1) 
obtained  interferograms  (Fig.  5)  which  implied  that  the  flow  was  not  planar. 

He  described  the  flow  disturbances  as  a  system  of  transverse  waves  in  the 
post-shock  relaxation  region  related  to  the  sinusoidally  curved  shock 
translational  front,  and  postulated  weak  forward-moving  oblique  shocks.  He 
observed  that  wave  geometry  was  not  consistent  in  identical  runs  although  these 
non-planarities  appeared  in  each  experiment.  It  was  found  that  no  references 
to  such  disturbances  were  made  by  other  researchers  investigating  similar 
ionized  argon  shocks  through  space  and  time-resolved  interferograms.  However, 
Bristow  directs  attention  to  Fig.  1  in  Ref.  23  and  a  system  of  weak  waves  in 
this  interferogra m  which  he  fails  to  identify  as  time-resolved.  This  moving- 
mirror  time- re solving  technique  gives  the  same  results  as  space-resolving 
techniques  only  when  the  flow  is  quasi- stationary,  that  is,  steady  in  the 
shock- fixed  coordinates,  so  that  time- independent  non-planarities  would  be 
faithfully  represented.  The  fact  that  all  perfect  gas  and  dissociating-gas 
flows  generated  by  this  facility  (Figs.  2  and  6)  were  always  planar  told 
him  that  this  was  an  ionization  phenomenon,  that  is,  the  electron  number 
density  ne  was  significant.  Bristow  suggested  that  the  much  larger  cross- 
section  of  this  facility  (18  x  10  cm  as  opposed  to  8.25  x  8.25  cm  in  Ref.  24 
and  5  x  5  cm  in  Ref.  23)  made  possible  an  ionization  resonance  effect  capable 
of  generating  such  transverse  waves. 

Morse  and  Ingard  (Ref.  7)  had  dealt  with  the  amplification  of  acoustic 
disturbances  in  an  electrically-heated  two- temperature  (Te  »  Ta)  plasma.  Since 
the  thermally-heated  ionizing  argon  plasma  under  consideration  here  has 
Ta  >  Te  in  the  relaxation  region  (Fig.  7a),  Bristow  conjectured  with  some 
experimental  but  no  theoretical  justification  that  this  temperature  difference 
could  also  cause  acoustic  amplification.  Measures  and  Belozerov  (Ref.  25) 
showed  that  for  strong  shocks  in  hydrogen  the  electron  and  heavy-particle 
temperatures  are  essentially  equal  in  relaxation  regions  (Fig.  7b).  Bristow 
reasoned  that  adding  some  hydrogen  to  the  argon  test  gas  would  cause  the 
difference  between  the  electron  and  heavy-particle  tenperatures  to  be 
significantly  reduced,  thereby  spoiling  the  acoustic  amplification  so  that 
the  flow  disturbances  would  not  appear.  The  addition  of  0.4$  of  H2  by 
pressure  to  the  argon  test  gas  in  fact  rendered  the  entire  flow  planar  (Fig.  4). 

It  was  found  that  ignoring  the  effect  of  the  hydrogen  inpurity  on  the 
equilibrium  fringe  shift  introduced  errors  smaller  than  those  from  other  sources, 
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and  the  investigation  continued  on  its  original  course.  Bristow  discussed 
several  sources  of  acoustic  disturbances  inside  the  shock  tube,  which, 
however,  should  give  rise  to  similar  non-planarities  in  perfect-gas  tod 
dissociating-gas  shocks  if  they  were  in  fact  the  cause.  The  only  plausible 
explanation  that  fits  the  experimental  observations  seemed  to  be  that  the 
postulated  acoustic-resonant  cavity  behind  the  shock  front  picked  up 
disturbances  from  sources  far  upstream  of  the  test  section  such  as  gas  inlet 
ports  and  amplified  them,  transmitting  them  down  the  shock  tube  into  the  tefct 
section  as  non-planarities  in  the  complex  reacting-gas  shock  wave.  In  the 
absence  of  a  resonant  cavity  disturbances  from  such  sources  die  out  rapidly, 
that  is,  as  x-^'2  where  x  is  the  distance  over  which  they  have  propagated 
(Ref.  2 6).  An  attempt  was  made  to  reduce  possible  effects  from  the  gas  ports 
to  the  rectangular  channel-,  but  it  met  with  little  success. 

The  addition  of  the  hydrogen  impurity  was  used  as  the  method  of  stabil¬ 
izing  the  shock-wave  flow  because  it  was  effective,  expedient  and  had  no 
detrimental  effects  on  the  flow  properties. 

Igra  (Ref.  8)  investigated  ionizing  argon  flow  over  an  expansion  corner 
by  causing  an  ionizing  shock  wave  to  pass  down  the  shock  tube  over  a  very 
long  (5.8m)  false  upper  wall  3.18  cm  thick  which  ended  at  the  test  section 
with  a  wedge  tappering  to  the  true  upper  wall.  Although  interested  in  the 
equilibrium  flow  behind  the  shock,  he  took  several  pictures  of  the  shock 
just  passing  over  the  corner.  They  all  showed  that  the  ionizing  argon  shock 
travelling  down  the  smaller  cross  section  was  always  planar.  Note  that  in 
this  instance  the  channel  height  was  no  longer  17.78  cm. 

Brimelow  (Ref.  9)  studied  shock-induced  boundary  layers  very  near  the 
shock-tube  bottom  wall  and  found  that  non-planarities  in  the  ionizing  argon 
shock  did  not  hinder  his  interferometric  measurements.  He  also  showed  that 
in  pure  argon  premature  ionization  behind  the  shock  front  '  occurred  close  to  the 
wall.  It  was  very  difficult  to  tell  if  this  phenomenon  was  present  when  the 
stabilizing  hydrogen  impurity  was  added.  No  link  was  made  between  this 
premature  ionization  near  (l  to  2  mm)  the  wall  and  the  shock-wave  non-planarities. 

Liu  (Ref.  11 )  conducted  an  analytical  study  on  the  effect  of  a  hydrogen 
impurity  on  the  shock  structure  of  ionizing  argon,  comparing  his  model- 
profiles  with  data  from  Ref.  9. 

Tang  (Ref.  10 )  investigated  experimentally  the  effects  of  different  types 
of  impurities  on  the  structure  and  the  non-planarities  of  ionizing  argon  shocks. 

He  found  that  oxygen,  helium,  and  sodium-chloride  had  no  observable  effect  on 
the  instabilities.  Only  the  addition  as  impurity  of  0.4£  H2  or  a  certain 
undetermined  amount  of  water  vapour  had  the  stabilizing  effect.  Tang  noted  that 
his  experiments  did  not  have  real-gas  state  variables  in  the  range  such  that 
(32p/3v2)s  <  0,  precluding  the  possibility  of  shock-wave  instability  associated 
with  the  violation  of  the  second  law  of  thermodynamics. 

In  bis  recommendation  for  future  study,  Bristow  suggested  that  an  analysis 
of  the  acoustic-wave  process  for  thermally-heated  plasma.  In  conjunction  with 
this,  it  was  considered  important  to  find  the  time  behaviour  of  the  non-plan&rities. 
The  fact  that  the  non-planarities  were  always  present  under  suitable  conditions 
although  the  wave  geometries  (e.g.,  the  shape  of  the  translational  shock  front) 
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were  not  reproducible  in  similar  runs  suggests  that  they  could  be  time 
dependent.  Normal  single-exposure  interferograms  gave  no  clue  as  to  the  motion 
of  the  sinusoidal-shaped  shock  front.  Also  it  was  not  clear  whether  these 
instabilities  occurred  across  the  width  of  the  shock  tube  as  well  as  the 
height,  that  is,  if  the  non-planarities  were  three-dimensional  or  only  two- 
dimensional. 

In  order  to  investigate  the  time  behaviour  of  the  non-planarities, 

G.  Doubilet*  used  a  framing  camera  to  obtain  five  interferograms  of  the  same 
shock  wave  in  ionizing  argon,  using  the  ruby  laser  in  the  normal-lasing  mode 
as  a  continuous  light  source.  Each  frame  took  0.5  microseconds  to  scan  and 
the  time  between  frames  was  adjustable;  one  photographic  plate  held  all  five 
12  mm  diameter  images.  In  contrast,  a  normal  single-exposure  interferometric 
image  of  the  test  section  is  about  10  cm  in  diameter.  The  laser  was  deficient 
as  a  continuous  light  source  because  in  normal  lasing  the  emitted  light 
consisted  of  pulses  about  0.3  microseconds  in  duration  and  random  in  time  and 
energy.  Figure  8a(i)  shows  a  trace  of  the  framing-camera  shutter  voltage  with 
time,  while  8a(ii)  shows  the  relative  intensity  of  light  emitted  from  the  ruby 
operated  in  the  normal-lasing  mode .  The  exposures  of  the  five  corresponding 
interferometer  frames  are  shown  in  8(b)  and  depend  on  the  strength  or  absence 
of  the  random  light  pulses  when  the  electronic  shutter  is  open.  The  two  most 
promising  experiments  are  shown  in  Figures  9a  and  9b  where  each  frame  has 
been  printed  separately  by  the  author  from  Doubilet' s  two  negatives  to  allow 
for  variations  of  exposure  of  the  same  type  as  seen  in  Figure  8b.  In  each 
frame,  the  translational  shock  front  can  be  seen  advancing  from  right  to  left, 
but  the  grain  caused  by  the  framing  camera  obscures  the  geometry  of  the  shock 
fronts.  Unfortunately ,  these  pictures  provide  no  information  about  the  time 
behaviour  of  the  non-planarities. 

Lapworth  (Ref.  26)  studied  the  evolution  of  shock-wave  perturbations  by 
placing  ahead  of  the  test  section  wedges  at  the  top  and  bottom  of  a  shock 
tube  and  taking  a  schlieren  picture  of  the  resulting  shock  wave  as  it  passed  the 
windows.  The  distance  between  the  wedges  and  the  windows  was  increased  between 
successive  pictures  to  observe  how  the  shock  evolved  with  distance  from  the 
wedges.  If  the  resulting  pictures  were  placed  side  by  side  according  to  the 
distance  from  the  wedges,  the  product  gave  the  illusion  of  being  one  picture 
consisting  of  many  exposures  of  the  same  shock  travelling  down  a  shock  tube  with 
transparent  sidewalls.  This  method  of  visualizing  time  behaviour  is  feasible 
only  if  both  the  shocks  and  the  observed  phenomenon  are  reproducible,  that  is, 
the  same  in  a  set  of  runs.  In  this  case  the  investigators  were  confident  that 
the  shock  strengths  and  the  subsequent  wave  geometries  were  consistent  over  a 
set  of  runs. 

The  above-mentioned  elegant  method  cannot  be  used  to  observe  the  time 
behaviour  of  the  present  phenomenon  for  two  reasons.  First,  the  source  of  the 
perturbations  is  unknown,  that  is,  the  researcher  has  little  or  no  control  over 
the  beginning  or  timing  of  the  disturbances  in  any  one  experiment.  Second,  the 
many  previous  interferograms  show  that  wave  geometry  (for  example,  the  shape  of 
the  translational  shock  front)  of  the  non-planarities  is  by.no  means  consistent 
between  similar  runs.  Thus,  the  method  employed  in  the  present  study  had  to 
ensure  that  multiple  exposures  were  obtained  of  the  same  shock  wave,  since  there 
was  no  known  reference  point  ahead  of  the  test  section  for  this  phenomenon . 
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Bristow  also  noticed  that  the  curved  shock  front  with  the  subsequent 
transverse  waves  appeared  similar  to  the  multi-headed  wavefronts  of  gaseous 
detonations  of  Hg, 02-mixtures,  as  reported  by  White  (Ref.  28).  The  structure 
and  evolution  with  time  of  these  wavefronts  are  well  understood  and  are 
illustrated  clearly  in  Fig.  5  of  Ref.  32.  Each  wavefront  is  a  shock  wave 
accross  which  runs  a  series  of  intersecting  triple  points,  each  wavehead 
consisting  of  a  length  of  the  shock  between  two  triple  points.  The  trajecto¬ 
ries  of  these  triple  points  visualized  by  allowing  the  passage  of  the 
pressure  pulses  and  slipstreams  associated  with  them  to  leave  imprints  on 
soot  deposited  on  the  shock-tube  wall.  It  was  suggested  that  this  method  be 
used  to  gain  some  understanding  in  a  similar  way  of  the  time  behaviour  of  the 
non-planarities  considered  in  this  study.  If  triple  points  or  some  other 
disturbances  existed,  leaving  transverse  traces  on  a  soot-covered  wall,  then 
a  cross-hatch  pattern,  such  as  that  Which  appears  in  Fig.  6  of  Ref.  32,  would 
be  expected. 

A  simple  relation  is  put  forward  in  Ref.  23,  connecting  the  average  size 
of  a  wavehead  with  the  detonation  velocity  and  the  reaction  induction  time. 

With  this  in  mind,  the  possibility  of  obtaining  an  empirical  formula  predicting 
the  number  of  wavelengths  in  the  sinusoidally  shaped  shock  front  was  investigated 

Therefore,  the  purpose  of  the  present  study  was  to 

1)  determine  the  time  behaviour  of  the  non-planar  shock  front, 

2)  determine  if  the  instabilities  are  also  three-dimensional,  and 

3)  see  what  can  be  learned  about  these  instabilities  by  the  soot-inprint 
technique. 

2.  EXPERIMENTAL  FACILITIES 


One  object  of  the  experimental  work  was  to  investigate  the  time-dependent 
behaviour  of  the  sinusoical  instabilities  of  the  translational  front  of  strong 
shock  waves  produced  in  argon  in  the  UTIAS  10  x  18  cm  combustion-driven 
hypervelocity  shock  tube.  This  was  accomplished  by  obtaining  in  one  picture 
severed  space-resolved  images  of  the  same  travelling  shock-wave  translational 
front  at  different  times.  Such  pictures  were  produced  by  a  multiple-exposure 
schlieren  system  based  upon  a  pulsed  ruby-laser  light  source  and  the  23  cm  dia 
field-of-view  Mach-Zehnder  interferometer.  This  instrument  is  installed  at  the 
shock  tube  test  section.  Many  modifications  (Ref.  1,19)  toward  the  improvement 
of  performance  and  operation  have  been  made  since  the  construction  (Ref.  22,26) 
of  the  facility. 

2.1  Shock  Tube 


The  UTIAS  hypervelocity  shock  tube  is  a  large  (15  m  long),  robust  facility 
built  to  simulate  high-enthalpy,  real-gas  flows,  generating  shock  speeds  in 
argon  of  up  to  8  km/s.  A  detailed  description  is  given  by  Boyer  (Ref.  22)  and 
an  excellent  update  is  found  in  Ref.  19.  Only  a  brief  description  will  be  given 
here,  including  modifications  done  by  the  author. 
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The  driver  section  consists  of  a  steel  cylinder  1.42  m  long  vith  an 
internal  diameter  of  15 «2U  cm  and  a  wall  10.l6-cm  thick.  Driver  pressures 
from  800  psi  to  6500  pse  were  obtained  by  filling  the  driver  vith  a 
combustible  mixture  of  7.5Jf  0g ,  20JS  H2  and  72.5/f  He  by  pressure  to  a  pressure 
p_£X,  and  igniting  this  mixture  by  impulsively  heating  a  0.38  mm  dia  tungsten 
wire  (strung  down  the  axis  of  the  driver)  vith  a  high-voltage  capacitor 
discharge.  The  combustion  products  sure  chiefly  hot  He  and  H^P  at  a  pressure  of 
about  6.7  times  Pmix*  constituting  a  near  idesLl  driver  gas  for  deflagrating 
shock-tube  operation  (Ref.  4).  Ten  percent  more  pressure  vas  added  to  ensure 
that  the  desired  driver  pressure  p^  vas  attained  and  the  diaphragm  did  burst. 
Thus  pm-[  -y  vas  related  to  the  anticipated  driver  pressure  p^  by 


=  (p^/6.7)  x  1.1 


(2.1) 


Each  diaphragm  vhich  vas  made  of  annealed,  type  SS  304  stainless  steel 
35.4  cm  in  dia.  Tvo  90°  grooves  vere  machine-scribed  on  the  driven  side  of  the 
diaphragm  at  right  angles  across  its  unsupported  diameter  to  ensure  petalling 


vhen  the  diaphragm  burst.  The  breaking  pressure  P^urst  of  t^ie  scribed 
diaphragm  vas  found  experimentally  by  Bristov  (Ref.  1)  to  be  given  by 

_  _  „ f4at)  fhln  . 

^burst  ~  1  d  J  [tj 


(2.2) 


vhere  t  -  total  diaphragm  thickness 

h  -  residual  thickness  in  scribed  grove 
d  -  unsupported  diameter  of  the  diaphragm  (=18.1  cm) 
a  -  ultimate  tensile  strength  of  SS  304  (=85,000  psi) 
n  -  fitted  constant  (=2.2  ±  0.1) 

’  1.01  vhen  t  a,  0.172"  (0.438  cm) 

and  K  -  ■  1.10  vhen  t  ^  0.108"  (0.295  cm) 

1.06  vhen  t  ^  0.062"  (0.158  cm) 

To  ensure  proper  petalling  of  the  steel  diaphragm,  the  condition  0.55  <  h/t  <  0.75 
vas  satisfied.  Repeatability  and  reliability  of  breaking  pressure  vas 
experienced  since  the  thicknesses  t  and  h  vere  held  vithin  rigid  machining 
tolerances.  To  achieve  a  given  shock  Mach  number  M  into  a  test  gas  at  pressure 
Pl  a  certain  driver  pressure  p^  is  needed.  Bri stoves  empirical  relation  (Ref.  l) 
betveen  p^,  p^,  and  Ms  for  argon  driven  by  the  combustion  process  outlined  above 
vas  used: 


ln(p. /p  )  =  0.422M  +  4.02  for  10  <  M  < 

4  1  S  S 


(2.3) 


For  example,  to  obtain  MB  =  16.15  in  2.50  torr  of  argon  a  driver  pressure 
of  p^  =  2452  psi  vas  required,  as  stipulated  by  Eq.  2.3.  A  diaphragm  of 
thickness  5  *  0.1805"  (0.458  cm)  vas  scribed  to  give  in  the  groove  a  residual 


— 'SfT' 


thickness  of  h  =  0.1015"  (0.257  cm),  so  that  it  would  burst  at  Pburst  *  2U52  psi 
(l68  atm)  =  pij,  per  Eq.  2.2.  To  produce  this  shock  wave,  p^  =  2.50  torr  of 
argon  was  admitted  intb  the  evacuated  channel,  and  into  the  driver  chamber  was 
admitted  30.2  (2.05  atm)  psi  of  Og;  He  gas  up  to  185  psi  (12.5  atm),  Hq  gas 
up  to  2 66  psi  (18.2  atm),  and  then  He  was  admitted  up  to  the  pressure  p^y  =  U03 
psi (27. 3  atm).  If  the  diaphragm  had  not  burst,  the  driver  pressure  would  have 
reached  a  maximum  of  about  6.7  x  Pmix  =  2700  psi  (l8U  atm)  after  combustion 
(Eq.  2.1).  As  it  was,  the  diaphragm  burst  at  about  2500  ±  50  psi  (170  i  3.U  atm), 
and  Mg  was  measured  at  the  test  section  as  16.17.  The  trace  of  the  driver 
pressure  with  time  is  shown  in  Figure  10a,  produced  by  a  pressure  transducer 
embedded  in  the  end  wall  of  the  driver,  opposite  the  diaphragm.  The  trace  starts 
when  the  capacitor  discharges  and  initiates  combustion,  and  shows  the  combustion 
pressure  rise  which  is  terminated  by  the  arrival  of  the  rarefaction  wave. 

Figure  10b  shows  the  pressure  record  in  the  test  section,  the  gauge  response 
time  (  10  us)  and  the  approach  to  the  equilibrium  flow  pressure  in  the  test 

section  followed  by  a  fairly  uniform  pressure  stall. 

2.1.2  Driven  Section 


The  15. 2U  m  long  driven  section  contains  a  0.457-m  circular  to  rectangular 
cross-section  transition  section,  a  1.22-m  test  section  located  13.56  m  from 
the  diaphragm  and  a  large  1  m3  shock-absorber  dump  tank.  The  channel  is  smooth 
on  all  four  sides  for  a  distance  of  11.6  m  ahead  of  the  20.32  cm  dia  test- 
section  windows,  9-53  cm  thick.  Upon  evacuation  of  the  driven  section  an 
ultimate  pressure  of  3  x  10“  torr  was  reached,  representing  an  impurity  of  200 
ppm  when  p1  =  1.5  torr  and  30  ppm  when  p^  =  10  torr.  Most  of  the  remaining  gas 
was  water  vapour.  The  argon  test  gas  had  a  quoted  impurity  level  of  less  than 
2  ppm.  The  pressure  p^  in  the  channel  was  measured  by  a  sensitive  oil 
manometer  giving  readings  from  0  to  &0  torr  by  the  formula 


_  1.069  +  (25  -  T)  x  9.5  x  10 


-1+ 


13.595 


x  H 


where  T  is  the  lab  temperature  in  cC,  H  is  the  oil  height  in  mm  and  p-^  is  in 
torr.  The  test  gas  temperature  T^  was  taken  as  the  reading  given  by  a  thermometer 
embedded  in  the  upper  wall  of  the  channel. 


The  arrival  of  the  shock- front  pressure  rise  was  detected  by  flushmounted 
piezoelectric  transducers  at  six  stations  called  D,  F,  G,  H,  I,  and  J  along  the 
channel  (Fig.  11a).  The  shock' a. arrival  at  D  started  five  counters  (No.  1  to 
No.  5)  connected  to  a  common  1  MHz  time  base,  and  the  arrival  of  the  shock 
front  at  each  of  the  remaining  five  stations  stopped  a  counter  in  turn.  The 
measured  travel  times  and  the  corresponding  distances  between  stations  were 
plotted  on  an  (x,t)  graph  with  x=0  at  D,  and  the  least-squares  second  order  fit 
through  these  point*  was  computed  giving  the  expected  travel  times  and  the 
shock  velocity  at  each  station  (Fig.  lib).  The  shock  Mach  number  quoted  for 
shock  waves  photographed  in  the  test  section,  and  used  in  subsequent  calculations 
for  that  experiment,  is  Ms  *  U/alt  where  U  is  the  expected  velocity  given  by 
the  second  order  fit  at  x  =  5*791  m  (station  M) ,  midway  between  I  and  J,  and 
ai  is  the  test  gas  sound  speed  given  by  a*  =  /yRTi ,  y  *  1.667  and  R  *  208.13 
a*/ s2-°K  for  argon. 
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The  average  deviations  between  times  observed  and  those  computed  by  the 
second  order  fit  (t  =  ax2  +  bx  +  c)  were  consistently  less  than  the  counter  reso¬ 
lution  of  ±1  us  over  the  whole  length  of  6.096  m  from  D  to  J  for  Mach  numbers 
from  9  to  17.  In  fact,  a  deviation  of  more  than  2  us  at  a  station  was  an 
indication  that  it  needed  servicing;  cleaning  and  recoating  the  piezoelectric 
transducer  with  a  thin  layer  of  silicone  rubber  usually  sufficed. 

As  a  check  an  additional  counter  (No.  7),  a  Hewlett-Packard  universal 
counter  (Type  5325A)  with  a  resolution  of  ±  0.1  us,  measured  the  time  taken  for 
the  shock  to  travel  from  I  to  J.  The  average  velocity  from  I  to  J  computed  from 
this  time  was  compared  with  the  computed  velocity  U  at  station  M  midway 
between  I  and  J  and  was  found  to  differ  by  only  a  few  percent  in  most  cases. 

Three  changes  were  made  in  this  study  affecting  the  accuracy  and  reli¬ 
ability  of  the  value  quoted  as  the  velocity  of  the  shock  wave  observed  in  the 
test  section. 

1)  The  three  Hewlett-Packard  (Type  373^A)  digital  counters  were 
assigned  to  stations  H,  I  and  J  which  were  closest  [2  ft  (0.6l  m),  1  ft 
(0.305  m)  and  1  ft,  respectively]  to  the  test  section  station  M,  while  the 
two  older  Racal  (Type  SA.1*5)  counters  with  observed  resolutions^ of  2  us 
and  4  us  were  assigned  to  stations  G  and  F,  farthest  from  M  [7  ft  (2.1U  m) 
and  H  ft  (3.36  m),  respectively].  Grouping  the  more  accurate  counters 
around  the  point  of  interest  was  expected  to  increase  the  reliability 

of  interpolations  at  M.  Previous  to  this  study,  the  least  accurate 
counters  were  assigned  to  I  and  J.  It  was  observed  that  thiB  change 
effected  a  significant  decrease  in  the  difference  between  the  velocity  at 
M  computed  from  the  fitted  curve  and  the  average  velocity  given  by  the 
time  measured  between  stations  I  and  J  on  counter  No.  7. 

2)  Previously,  the  order  of  the  fitted  curve  was  the  lowest  for  which 
the  RMS  deviation  between  observed  and  expected  times  was  less  than  2  us. 
Thus  an  incorrect  time  caused  by  a  malfunctioning  pressure  transducer 
could  lead  to  high-order  fitted  curves  that  gave  nonphysical  results. 
Locating  and  eliminating  the  erroneous  datum  required  a  manual  work- 
through  of  the  data.  This  disadvantage  waB  aggravated  by  additional 
errors  caused  by  assigning  to  the  stations  I  and  J  the  old  Racal  counters 
( see  .previous  paragraph)  which,  when  set  to  record  the  same  time, 
typically  disagreed  by  U  us  or  more,  even  when  using  the  same  time  base. 

For  this  study,  a  second  order  least-squares  curve  was  fitted 
through  the  data  points  and  the  shock  velocity  quoted  was  the  expected 
velocity  computed  from  this  curve  at  M.  The  choice  of  a  second  order 
curve  seemed  physically  appropriate  because  the  expected  shock  velocity 
would  decrease  monotonically,  reflecting  the  attenuation  of  the  shock  by 
viscous  forces.  This  was  born  cnxt  by  the  fact  that  the  data  consistently 
fit  the  second  order  curve  very  well  for  Ms  from  3.0  to  18.3.  The  computer 
program  which  did  the  calculation  included  the  deviations  of  all  the 
observed  times  from  the  fitted  times  so  that  an  errant  datum  could  be 
easily  spotted  and  quickly  removed.  These  deviations  did  not  exhibit 
any  systematic  non-random  pattern  in  this  Mach  number  range,  indicating 
further  the  adequacy  of  the  second  order  fit. 
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3)  In  Ref 8.  1  and  10  the  shock  velocity  quoted  was  the  average 
velocity  calculated  from  the  time  interval  measured  hy  counter  Ho.  7, 
and  this  was  check  against  the  expected  velocity  at  M;  however,  this 
method  depended  upon  the  good  working  order  of  the  pressure  transducers 
at  I  and  J  and  upon  the  shock  wave  being  quasi-stat ionary.  Reference  19 
refers  to  a  method  of  obtaining  a  value  for  U  that  depends  upon  the 
average  velocity  between  H  and  X.  After  observing  the  magnitude  of  the 
non-planarities  of  some  of  the  shock  fronts  obtained  before  (Refs.  1,10), 
it  was  felt  that  possible  shock-front  oscillations  would  render  unreliable 
the  average  velocities  obtained  from  the  time  measurements  along  only  a 
few  feet.  Thus  the  second  order  fitting  curve  was  used  to  smooth  out 
any  potential  variations  in  observed  times  due  to  shock-front  oscillations. 
Using  the  expected  value  of  U  at  M  would  also  reduce  the  dependence  of  U 
on  the  working-condition  of  a  few  key  counters  and  pressure  transducers. 

It  must  be  used,  of  course,  when  a  model  is  installed  in  the  test  section, 
since  the  shock  may  be  changed  by  the  time  it  reaches  station  J. 

2.1.3  Interferometer  and  Light  Source 

The  diagnostics  of  the  gas  flow  in  the  UTIAS  hypervelocity  shock  tube  are 
performed  with  the  aid  of  a  23-cm  dia  aperture  Mach-Zehnder  interferometer. 
Details  of  the  design  and  operation  of  thiB  device  are  found  in  Ref.  20. 

Up-dated  alignment  procedures  are  given  in  Ref.  19.  A  short  description  of  the 
optical  components  involved  will  facilitate  a  discussion  of  the  application  of 
interferometry  to  the  present  study. 

The  interferometer  whose  dimensions  are  given  in  Fig.  12  is  arranged  such 
that  the  upper  arm  contains  the  test  section  which  is  aligned  properly  with  the 
axis  of  the  device  by  a  continuous  He-Ne  gas  laser.  The  lower  arm  contains  a 
compensating  chamber  with  identical  windows  and  test  gas.  The  peripheral  optics 
are  shown  schematically  in  Figure  13.  The  light  source  for  this  instrument  is 
a  TRG,  Model  10U,  pulsed  ruby-laser  system.  The  8  mm  dia  beam  of  parallel  light 
is  focussed  by  a  lens  to  the  desired  source  configuration  at  f^  (point  for 
interferometry  and  shadow,  narrow  rectangle  for  schlieren).  The  light  which 
passes  through  the  test  section  as  a  parallel  test  beam,  due  to  defocussing  by 
parabolic  mirror  1  and  reflection  by  splitter  1  and  mirror  1,  is  focussed  by 
parabolic  mirror  2  to  f2  and  passes  into  the  camera.  For  interferograms  the 
compensating  beam  consisting  of  that  light  passing  through  splitter  1  is 
reflected  by  mirror  2  and  splitter  2,  interfering  with  the  test  beam.  For  .. 
shadow  and  schlieren  pictures  of  events  inside  the  test  section,  the 
compensating  beam  is  blocked.  The  test  section  object  plane  is  imaged  at  the 
camera' 8  ground  glass  by  parabolic  mirror  2  and  the  camera  lens  with  a  net 
magnification  of  about  0.5. 

A  short-duration  pulse  of  light  is  obtained  from  the  ruby  laser  by 
0-switching  the  laser  cavity  with  a  Pockels  cell  electro-optical  shutter  and 
polarizer.  A  retarding  voltage  of  3  kV  is  applied  to  the  Pockels  cell 
preventing  lasing  until  the  picture  is  desired.  At  this  time,  the  foltage 
is  removed,  that  is,  the  laser  is  Q-switched,  and  lasing  converts  the  stored 
energy  intg  a  30-ns  pulse  of  coherent,  monochromatic  light  of  wavelength 

«  69 U3  A.  A  second-harmonic  generator  (SHG)  provides  a  shorter  wavelength 
X2  *  31*71.5  A.  The  total  optical  energy  output  is  0.8  Joules . 
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The  distance  between  the  desired  shock  position  and  station  I  divided 
by  the  predicted  (see  Eq.  2.3)  shock  speed  is  used  as  the  time  delay  betveen 
the  detected  arrival  of  the  shock  at  I  and  the  triggering  of  the  Q- switch. 

The  difference  between  the  observed  and  desired  shock  position  depends  on  the 
predictibility  of  the  shock  speed,  and  is  usually  within  a  few  centimeters 
over  the  23-cm  field  of  view. 

The  flash-lamp  pumping  time  tp  is  the  time  between  the  triggering  of  the 
laser  flash  lamp  and  the  Q-switching,  and  normally  is  designed  to  be  900  us 
for  optimal  single-exposure  photos.  A  different  value  of  tp  ■  500  us  was 
used  for  some  of  the  present  experiments.  The  station  which  is  closest  to 
the  test  section  and  which  is  expected  to  detect  the  shock  at  a  time  at  least 
tp  before  the  time  that  the  shock  is  calculated  to  reach  the  desired  position, 
initiates  a  time  delay  unit  which  accounts  for  the  excefcs  time  over  t_  and 
starts  the  flash-lamp  pumping.  Counter  No.  6  records  the  actual  time  interval 
tp  between  the  flash-lamp  trigger  and  the  Q-switch  which  is  determined  mainly 
by  the  actual  shock  speed.  The  measured  t_  for  an  expected  900  us  ranged 
from  850  to  950  us  with  no  consequential  loss  in  picture  quality.  The 
measured  tp  for  an  expected  500  us  ranged  from  U60  to  5^0  us  and  had  an 
effect  on  the  strength  of  the  multiple  pulses.  This  showed  up  in  the  number 
of  visible  translational  shock-front  images  in  one  multiple-exposure 
schlieren  picture. 

Focussing  of  the  test  section  object  plane  is  facilitated  by  providing 
cross-wires  on  the  windows  which  can  be  sharply  defined  on  the  ground  glass 
and  film  plates.  Where  the  wires  cross  can  be  used  as  absolute  reference 
points  whose  separation  can  be  measured  '  .25  mm)  so  that  an  accurate  scale 
for  observed  phenomena  can  be  provided  in  the  final  pictures.  The  wires  can 
also  be  strung  in  the  x  and  y  directions  (correct  orientation  ±0.01  rad)  to 
conveniently  establish  the  channel  axes  on  the  final  picture  (see  Fig.  l). 

In  the  camera  the  light  passes  through  the  camera  lens  and  through  a 
Xj-line  filter,  which  effectively  filters  out  the  second  harmonic  and  the 
broadband  radiation  from  the  ionizing  argon  shock,  and  onto  the  film  plate 
marked  X^.  When  dual-wavelength  interferograms  are  desired,  part  of  the  light 
is  reflected  by  a  beam  splitter  and  is  folded  by  a  mirror  down  the  second  arm 
of  the  camera,  and  passes  through  a  Xg-line  filter  to  expose  the  film  plate 
marked  Xg.  Kodak,  Royal  Pan  X  (1250  ASA),  10  cm  x  13  cm  plate  film  was  used. 
Upon  exposure  of  the  plates,  which  included  a  no-flow  set  taken  immediately 
before  each  experiment,  they  were  processed  in  Kodak,  DK  50  developer  for  7 
minutes,  stop  solution  for  15  seconds,  and  then  fixer  for  20  minutes. 

2.1.U  Accuracy  of  Measurements 

The  initial  temperature  is  measured  directly  from  a  standard  mercury- 
bulb  thermometer  with  an  uncertainty  of  ±0.1°K.  Thus  the  absolute  error 
associated  with  T^  is  E(T^)  a  0.2°K. 

The  initial  pressure  is  measured  indirectly  from  a  sensitive  oil 
manometer  from  which  the  oil  height  is  observed.  The  absolute  error  asso¬ 
ciated  with  any  initial  pressure  from  0.2  to  28  torr  is  less  than  0.1  torr. 


The  incident  shock-wave  Mach  Humber  Mg  must  be  calculated  from  the 
shock  velocity  and  initial  sound  speed  in  the  argon  gas.  The  shock  speed  is 
obtained  from  a  measured  distance  S,  traversed  in  a  measured  time  t,  while  the 
sound  speed  is  calculated  from  the  initial  temperature  T^  mentioned  above. 

The  uncertainties  in  the  distance,  time  interval  and  temperature  measurements 
cause  a  maximum  relative  error  associated  with  M  which  can  be  written  as 
(Ref.  19) 


E(M  ) 

-t; - -  (1.05  x  M  +  10.18)  x  10”  for  argon. 

MS 

s 

At  M  *  10,  the  relative  error  is  2 1  and  E(M8)  *  0.2  whereas  at  MB  *  18,  the 
relative  error  is  2.9Jt  and  E(Mg)  =  0.52.  The  error  E(Mg)  is  the  maximum 
possible  error  and  in  most  cases  vould  actually  be  much  less. 

As  discussed  in  the  previous  section,  the  method  of  calculating  the 
average  shock  velocity  between  I  and  J  by  measuring  the  time  interval  the 
shock  takes  to  traverse  the  short  distance  (S=6l  cm),  was  not  used  in  this 
study.  This  was  Justified  also  by  the  fact  that  the  reading  of  the  interval  t 
on  counter  No.  7  was  often  spoiled  by  circuit  noise  due  to  the  large  number  of 
free  electrons  (high  n*)  inside  the  steel  tube.  In  such  circumstances  it  was 
in  fact  necessary  to  use  the  expected  velocity  at  M  provided  by  the  second 
order  fitted  curve  for  the  shock  velocity.  When  a  good  reading  for  t  was 
obtained,  the  shock  velocities  derived  by  these  two  methods  differed  by  lesB 
than  1  /(,  the  latter  always  giving  a  greater  velocity. 

Bristow  (Ref.  l)  states  that  the  shock  oscillations  must  be  an  ionization 
phenomenon,  appearing  when  the  maximum  electron  number  density  ne  E  becomes 
significant.  The  temperatures  Ta  and  Te  also  seem  to  be  involved!  To  discover 
what  range  of  values  of  ne  „  and  temperature  T&  (reflected  in  MB)  were 
available  in  this  facility!  the  equations  in  Sect.  2.1  and  a  thermodynamic 
routine  were  used  to  calculate  Mfl  and  ne  E  for  all  the  values  of  p^  and  p> 
practical  in  this  shock  tube.  The  result  iB  plotted  in  Fig.  14  for  T^  *  298°K. 
The  driver  pressure  is  limited  to  p^  *  6500  psi  for  safety.  Large  values  of 
p^  and  ph  give  post-shock  pressures  too  great  for  the  test-section  windows 
(p2  *  100  psi).  Very  low  values  of  p^  give  densities  too  small  to  trigger  the 
shock  detectors.  Note  that  a  maximum  value  of  ne  E  is  obtained  for  a  given 
Pi,  for  a  Mach  number  of  about  1U.5.  An  absolute  maximum  value  for  this  shock 
tube  operation  of  ne  E  *  2.U8  x  lO1^  cm”3  is  possible. 

2.2  Multiple-Pulse  Light  Source 

2.2.1  Pulsed  Laser  in  Combined  Mode 


A  multiple-exposure  schlieren  system  required  a  multiple-pulse  light 
source  which  was  realized  in  a  ruby  laser  operated  in  a  combined  giant-pulse/ 
normal-lasing  mode  For  a  full  description  of  the  ruby  laser  and  discussions  on 
modes  of  operation,  see  Ref.  29.  Here  a  simple  description  will  reveal  the 
behaviour  behind  the  combined  mode  used  specifically  for  this  application. 

The  flash  lamp  inside  the  laser  irradiates  the  ruby  crystal  rod  for  about 
1.5  ms  and  the  absorbed  radiation  causes  population  inversion.  Normal  lasing, 
in  the  ruby  starts  as  soon  as  the  inversion  level  permits  it,  and  the  inverted  level 
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is  depleted  very  quickly  {fraction  of  a  ys),  giving  a  small  pulse  of  coherent 
monochromatic  light.  Consequently,  vhile  the  flash  lamp  is  on,  small  pulses 
random  in  time  and  energy  are  emitted  (Fig.  8a).  In  the  giant-pulse  mode,  the 
lasing-cavity  resonance ,  or  Q,  is  spoiled  hy  directing  incipient  radiation  away 
from  the  end  mirrors  until  the  desired  instant  when  the  Q-spoiling  is  switched 
off.  During  the  time  interval  tp  between  the  onset  of  discharge  through  the 
flash  lamps  and  the  Q-switching,  the  inversion  level  in  the  crystal  has  been 
allowed  to  build  up  to  a  level  much  larger  than  that  obtained  during  normal 
lasing.  After  the  Q-spoiling  is  removed,  lasing  occurs  until  the  inversion 
population  is  depleted,  dumping  in  about  30  ns  the  stored  energy  into  a 
"giant”  pulse  of  coherent  monochromatic  light. 

If  the  flash  lamp  is  still  on  after  Q-switching,  i.e.,  tp  <  1.5  ms, 
then  normal  lasing  occurs  after  the  giant  pulse.  The  energy  in  the  giant 
pulse  and  the  total  energy  in  the  normal  lasing  pulses  depend  upon  tp. 

Setting  tp  =  0  ys  disables  the  Q-switch  and  gives  normal  lasing  with  no  giant 
pulse.  Setting  tp  =  900  ys  gives  the  largest  giant  pulse  with  insignificant 
normal  lasing;  this  is  the  setting  used  to  take  interferograms  at  this  facility, 
in  order  to  get  one  set  of  bright  fringes  with  no  "ghost"  fringes  (Ref.  19). 

For  this  study  it  was  found  that  operating  the  Q-switch  at  an  intermediate 
value  of  t  produced  a  moderate-strangth  giant  pulse  and  several  very-short- 
duration  normal-lasing  pulses  of  usable  energy.  The  optimum  value  was 
tp  =  500  ys,  and  was  used  in  each  of  the  multiple-exposure  schlieren  experiments. 
In  practise  tp  varied  from  U6o  ys  to  5^0  ys;  the  reason  for  this  iB  given  in 
Sect.  2.1.2. 

The  reason  for  the  existence  of  an  optimum  value  of  t  is  clear  if  these 
simplified  arguments  are  considered.  Inductances  in  the  discharge  circuit  and 
a  finite  heating  time  in  the  flash  lamp  cause  the  maximum  power  output  of  the 
lamp  to  occur  about  500  ys  after  the  discharge  is  initiated,  that  is,  at  tp  *  500 
ys.  The  energy  in  the  giant  pulse  depends  upon  the  accumulated  energy  which  is 
reflected  in  the  area  under  the  flash-lamp  power  output  curve  corresponding  to  the 
pumping  period.  However,  once  the  giant  pulse  empties  the  energy  stored  in  the 
ruby  crystal,  the  energy  in  the  normal  lasing  pulses,  depends  upon  the  energy 
transfer  rate,  i.e.,  the  height  of  the  curve.  Hence  optimum  multiple-pulse 
output  is  achieved  when  energy  transfer  to  the  ruby  crystal  is  greatest,  about 
500  ys  after  discharge  initiation. 

For  shock  speeds  of  about  5  km/s  the  shock  wave  is  in  the  field  of  view  for 
about  50  ys.  The  light  output  of  the  laser  for  the  50  ys  following  a  giant  pulse 
with  tp  =  500  ys  was  obtained  by  measuring  on  an  oscilloscope  the  current  through 
a  photodiode  exposed  to  the  light.  Figure  15  shows  such  an  oscilloscope  trace. 

The  giant  pulse  p-^  is  an  order  of  magnitude  stronger  than  the  other  pulses,  and 
P2  is  the  second  strongest.  As  shown  by  this ‘and  all  other  traces  as  well  as 
all  schlieren  pictures,  P2  follows  P-j_  by  the  same  time  interval — about  7-3  ys. 

An  explanation  for  this  probably  involves  details  of  the  Pockels-cell  Q-switching 
process.  The  remaining  pulses  were  random  in  energy  and  quite  unreliable ,  in 
contrast  to  P2<  Over  about  1*0  experimental  runs  P-^  produced  a  clear  image  each 
time,  P2  produced  a  usable  image  in  nine  tenths  of  them,  and  a  third  image  caused 
byPg  was  seen  in  less  than  one  third  of  them. 

Unlike  the  normal  pulses  shown  in  Figure  8a  where  tp  *  0  ys,  those 
immediately  following  the  giant  pulse  shown  in  Fig.  15  had  very  short  durations. 
This  observation  was  reinforced  by  the  Bharpness  of  the  schlieren  images  caused 


•  m-.,, 


by  these  light  pulses.  At  this  shock  speed  no  perceptible  blurring  of  the  shock 
front  image  implied  that  the  pulse  durations  were  less  than  0.1  us . 

The  time  interval  between  P^  and  Pg  was  not  adjustable,  but  was  sufficient 
for  the  present  purpose.  In  fact,  it  was  very  convenient  for  two  reasons.  First, 
this  time  interval  was  long  enough  that  the  images  of  the  shock  wave  produced  by 
them  were  definitely  different,  that  is,  a  change  in  shape  was  clearly  seen. 
Secondly,  the  interval  was  short  enough  that  the  relationship  between  the  two 
images  was  clear.  The  presence  of  a  third  or  even  fourth  image  tended  to 
confirm  this. 

Thus  a  multiple-pulse  light  source  capable  of  determining  with  a  certain 
degree  of  confidence  the  time  behaviour  of  the  instabilities  was  realized 
simply  and  cheaply  merely  by  changing  the  operational  parameter  tp  using  the 
giant-pulse  mode  of  the  existing  ruby  laser.  Also,  it  was  fortunate  that  the 
light  pulses  used  were  monochromatic  of  wavelength  =  69U3  A  so  that  a  X^ 
notch  filter  could  pass  the  light  from  the  light  source  but  essentially  reject 
the  radiation  emitted  by  the  hot  plasma  behind  the  strong  shock  wave, 
before  it  reached  the  camera. 

2.2.2  Laser/Spark  Light  Source 

At  the  beginning  of  the  experimental  work  it  was  felt  that  two  separate 
controlled  light  pulses  could  be  obtained  from  the  ruby  laser  and  a  simple 
spark  light  source,  each  contributing  one  pulse.  The  time  between  pulses  was 
to  be  variable  to  provide  two  images  of  the  shock-wave  translational  front 
separated  by  an  adjustable  time  and,  hence,  distance.  The  incumbent  laser 
was  to  remain  unchanged*  producing  a: Cl 8  mm  dia  beam  of  parallel  light  con¬ 
verged  by  a  lens  to  a  point  at  f-j_  to  act  as  a  point  source  as  required  by 
the  interferometer  (Ref.  20) (Fig.  13).  The  spark  light  source  was  designed 
so  that  the  electrodes  straddled  the  point  f-j_,  producing  a  good  approximation 
to  a  point  source  at  fj  while  allowing  the  laser  light  to  pass  through  in  its 
turn.  The  spark  source  designed,  built  and  tested,  by  A.  Morte  had  a  0.01  yf 
capacitor  and  was  operated  at  about  15  kV,  thus  stored  about  1.0J,  of  which 
only  a  small  fraction  was  emitted  as  broadband  light  (Ref.  17).  Likely  less 
than  10“ 3j  of  this  was  collected  and  used  to  produce  a  photographic  image.  In 
contrast,  the  laser  power  supply  stored  1  kJ.  Of  that,  about  10  J  was  emitted 
as  radiation  by  the  flash  lamp,  and  about  1  J  was  emitted  by  the  ruby  crystal 
as  a  30-ns  pulse  of  coherent  monochromatic  light  of  wavelength  X^  with  a 
bandwidth  of  about  0.1  A  (Ref.  29). 

The  first  system  consisted  of  the  usual  Pockels  cell  triggering  system 
for  single-exposure  interferograms ,  to  which  was  added  only  a  time  delay  unit 
adjustable  in  the  range  0  to  100  ys  and  the  triggerable  spark  light  source.  The 
delay  unit  was  started  when  the  Pockels  cell  was  triggered,  that  is,  when  the 
laser  pulse  occurred,  and  the  delay  unit  triggered  the  spark  source  after  the 
set  time  delay. 

The  light  output  of  the  spark  was  Just  great  enough  to  give  a  good 
shadowgraph  image  on  a  Kodak  Royal  Pan  X  negative  plate  of  non-radiating  shock 
waves  in  the  absence  of  the  X,  notch  filter  which  is  a  necessity  for  obtaining 
interferograms  of  radiating  shocks  with  the  laser.  The  spark  emitted  a  light 
pulse  with  a  duration  less  than  50  ns  and  gave  sharp  images  of  shock  waves  up 
to  about  Ms  »  7  in  argon,  i.e.,  shock  speed  of  about  2.5  mm/ys.  However, 
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above  Mg  =  7  significant  radiation  from  the  plasma  behind  the  shock  overpowered 
the  spark  light.  Since  both  sources  were  emmitting  broadband  light  and  were  of 
comparable  brightness,  the  spark  image  could  not  be  seen  nor  separated,  and  the 
picture  was  spoiled. 

The  author  feels  that  this  two-exposure  system  would  be  very  useful  in 
studying  many  other  problems  based  on  evolving  non-radiating  shock  waves  in  the 
UTIAS  hypervelocity  shock  tube,  since  it  was  easy  to  set  up  and  reliable  in 
operation.  The  spark  source  was  subsequently  not  used  in  the  present  study 
only  because  the  phenomenon  considered  involved  radiating  shocks,  to  which  the 
above  mentioned  system  was  therefore  not  applicable.  Installing  a  capacitor  of 
higher  value  would  certainly  increase  the  spark  energy  and  hence  the  usefulness 
of  this  simple  yet  effective  auxiliary  spark  light  source. 

2.2.3  Discussion  of  Other  Light  Sources 

Reference  30  describes  the  design  and  operation  of  an  eight-channel  high¬ 
speed  multiple- spark  camera  which  gave  eight  separate  schlieren  pictures 
arranged  in  a  circle  on  a  U"  x*  5”  plate.  However,  the  energy  of  each  spark 
was  about  the  same  as  that  of  the  spark  source  described  in  Sec.  2.2.2,  and 
for  the  same  reason  the  light  source  could  not  be  used  here.  This  camera  could 
produce  eight  good  quality  schlieren  photos  separated  by  adjustable  times,  at 
the  cost  of  a  complex  triggering  system.  If  only  two  exposures  of  shock  waves 
will  suffice,  the  pulsed-laser/spark-source  system  has  the  great  advantage  of 
economy,  convenience  and  reliability. 

A  previous  spark  source  (Ref.  20)  used  alone  with  this  interferometer 
employed  a  very  high-energy  spark  with  a  notch  filter  to  obtain  a  strong  pulse 
of  light  with  a  very  narrow  frequency  distribution,  similar  to  that  of  a  laser 
beam  but  not  as  narrow.  A  similar  filter  was  placed  in  the  beam's  path  after 
the  test  section  to  filter  out  most  of  the  broadband  shock  radiation,  letting 
through  the  source  light  which  contained  more  energy  than  the  filtered  shock 
radiation.  Because  the  capacitance  of  the  discharge  circuit  was  so  high 
(22  yf),  the  pulse  duration  was  very  long  (5  ys)  and  the  pulse  had  to  be 
shuttered  to  less  than  0.1  ys  with  a  Kerr  cell.  This  system  was  later 
replaced  by  the  ruby  laser  and  the  Xj  notch  filter. 

Neither  the  high-power-spark  source  nor  the  Ker  cell  were  available,  so 
the  simple  spark  source  discussed  previously  was  tried.  Unfortunately,  none 
of  the  character! sties  of  the  instabilities  appeared  at  Mach  numbers  at  which 
the  spark  gave  a  usable  picture,  i.e.,  below  Mach  7-  Hence  the  simple  spark 
light  source  could  not  be  employed. 

The  ideal  additional  light  source  would  be,  of  course,  another  ruby  laser 
capable  of  giant-pulse  operation,  which  would  make  feasible  two-exposure 
photography  of  radiating-gas  flows.  Unlike  the  first  two  light  pulses  Pj_  and 
P2  obtained  from  the  ruby  laser  in  the  combined  giant-pulse  normal-lasing  mode 
discussed  above,  two  lasers  would  provide  light  pulses  of  adjustable  time 
separation  and  comparable  strengths.  This  would  obviate  the  need  for  the 
special  modification  of  the  schlieren  method  discussed  in  Sec.  2.3.2.  However, 
all  thie  ruby  lasers  available  on  loan  had  for  the  giant-pulse  Q-switch  a  rotating 
mirror  which  has  an  uncontrollable  time  delay  between  the  Q-switch  trigger  and 
laser  light  output  of  anywhere  from  0  to  125  ys.  These  lasers  were  therefore 
of  little  practical  use,  because  in  this  situation  where  the  shock  travels  at 
up  to  U.5  mm/ys  there  should  be  a  delay  between  trigger  and  light  pulse  of  no 
more  than  tenths  of  a  microsecond,  as  is  the  case  with  the  Pockels  cell.  Hence, 
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the  existing  laser  was  used  for  generating  two  pulses,  although  the  separation 
time  was  fixed  at  7  us.  As  pointed  out  previously,  this  compromise  presented  no 
serious  drawback  to  the  achievement  of  the  purpose  of  the  present  work. 

2.3  Optical  Techniques 

This  section  describes  how  the  schlieren  method  was  applied  to  the  space- 
resolved  multiple-exposure  visualization  of  shock  waves.  It  tells  why 
schlieren  was  chosen,  how  schlieren  was  adapted  to  multiple-exposure  photo¬ 
graphy,  and  how  an  optimal  picture  was  obtained. 

2.3.1  Space  vs  Time-Resolved  Photography 

The  non-planar  translational  front  of  strong  shocks  in  argon  was  first 
observed  in  interferograms  obtained  with  the  UTIAS  Mach-Zehnder  interferometer 
(Refs.  2,9).  The  theory  and  application  of  the  interferometer  are  discussed  in 
Refs.  l6  and  15,  but  a  few  points  pertaining  to  this  study  are  emphasized  here. 

The  interferometer  gives  a  circular  side  view  of  the  shock  tube  channel  which 
is  17.78  cm  (7.00")  high,  through  two  23  cm  (8")  diameter  9.515  cm  (3.75") 
thick  windows,  thus,  a  two-dimensional  space-resolved  image  of  the  test  section 
is  obtained,  as  opposed  to  a  time-resolved  image  produced  by  a  rotating-mirror 
camera  by  passing  the  image  of  a  narrow  slit  over  the  film.  Time-resolved 
pictures  are  most  meaningful  when  the  phenomenon  in  the  shock  tube  is  quasi¬ 
stationary  and  one-dimensional.  Space-resolved  interferograms  can  confirm 
one-dimensionality  by  showing  that  the  fringe  shifts  and  hence  flow  properties 
are  independent  of  y,  i.e.,  the  height.  In  the  present  case  the  shock  is 
two-dimensional,  in  the  x  and  y  directions,  so  the  interferometer  which 
resolves  the  x  and  y  directions  is  the  right  tool  to  use  to  observe  these  non- 
planar  effects.  Dependence  of  flow  variables  on  z  is  small  and  does  not 
significantly  effect  the  pictures  except  that  sometimes  the  translational  front 
has  the  appearance  of  a  twisted  ribbon,  as  seen  in  Fig.  3. 

To  observe  the  time  dependence  of  the  shock  front  the  two-dimensional  image 
of  the  shock  front  must  be  resolved  in  time.  This  cannot  be  done  continuously 
on  a  single  photographic  plate.  However,  since  the  shock  wave  is  travelling, 
discrete  pulses  of  light  through  the  test  section  will  cause  images  of  the 
translational  shock  front  to  appear  at  different  places  on  the  film  plate. 

Thus  a  multiple-exposure  picture  would  show  some  evolution  of  the  non-planarities 
of  the  shock  front.  The  multiple  pulse  provided  by  the  ruby  laser,  as  described 
in  Sec.  2.2.1,  was  used  for  this  purpose. 

2.3.2  Interferometry  vs  Schlieren  Method 

Several  multiple-exposure  interferograms  were  taken  to  test  the  technique, 
but  the  results  were  dissappointlng.  Giant  pulse  filled  the  entire  field 
of  view  with  bright  fringes,  and  the  additional  faint  fringes  produced  by  pg 
were  very  difficult  if  not  impossible  to  detect.  Figure  l6  is  a  multiple- 
exposure  interferogram  of  a  shock  wave  travelling  from  right  to  left  in  argon 
with  Ms  =  lU,3,  p^  =  10.2  torr,  T^  =  298.56°K,  and  tp  =  494  us.  The  fringes 
produced  by  the  pulse  are  clear,  and  the  shock-front  shape  is  evident  at  A. 

The  shock  front  for  the  second  exposure  should  appear  near  B,  but  the  P2  fringes 
are  barely  discernible,  giving  no  useful  information  about  the  time  behaviour  of 
the  phenomenon. 
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Multiple-exposure  interferograms  would  be  more  feasible  if  the  light  pulses 
were  of  lesser  and  comparable  energies  as  would  be  the  case  with  two  ruby- 
lasers.  Since  another  ruby  laser  with  a  Pockels  cell  was  not  available, 
it  was  concluded  that,  even  though  interferometry  was  the  most  useful  single¬ 
exposure  space-resolving  visualization  method,  it  could  not  be  used  with  the 
present  system. 

Since  only  the  translational  shock  front  was  initially  of  interest,  the 
schlieren  method  (Ref.  21)  of  visualizing  density  gradients,  which  would  show 
primarily  the  very  high  density  gradient  of  the  translational  shock  front,  was 
the  best  method  to  use  in  this  case.  Thus  it  was  expected  that  each  pulse  of 
light  which  passed  through  the  test  section  when  the  shock  was  present 
would  leave  on  the  photographic  plate  the  ambient  illumination  and  the 
schlieren  image  of  only  the  translational  front  of  the  shock.  The  image  of 
the  translational  front  would  appear  as  a  thin  vertical  line,  straight  if  the 
shock  front  was  planar,  or  curved  or  kinked  if  otherwise.  Hence,  when  several 
distinct  pulses  of  light  passed  through  the  test  section  while  the  shock  was 
travelling  in  it,  corresponding  images  of  the  translational  front  would  fall  on 
the  photographic  plate  in  different  places,  showing  the  position  and  geometry 
of  the  shock  front  at  the  times  at  which  the  light  pulses  occurred.  Each  thin 
vertical  image  would  be  displaced  horizontally  from  the  previous  one  by  a 
distance  determined  by  the  time  between  light  pulses  from  the  laser  and  the 
velocity  of  the  shock  wave.  Of  these  three  quantities,  the  shock  velocity  and 
image  displacement  were  accurately  measurable.  The  time  between  B-j_  and  Bo 
calculated  in  this  manner  could  be  compared  with  those  observed  on  the  traces 
of  laser  light  output  (Fig.  15).  This  time  was  consistently  observed  and 
calculated  to  be  7.3  i  0.5  ys.  In  addition,  it  was  found  that  the  images  of 
flow  phenomena  besides  the  shock  front  did  appear,  but  assigning  them  to  the 
correct  shock  front  image  was  facilitated  by  referring  to  single-exposure 
interferograms  and  schlieren  pictures  of  similar  experiments. 

2.3.3  Multiple-Exposure-Schlieren  Technique 

The  interferometer  is  readily  available  for  employing  the  schlieren  method 
to  visualize  density  gradients  in  the  shock  tube,  requiring  only  slight 
modifications  to  the  peripheral  optics  to  obtain  optimal  results.  Figure  12 
shows  the  particular  arrangement  used  in  this  study  (see  Sec  2.1.3). 

To  obtain  dual-wavelength  schlieren  pictures,  Igra  (Ref.  8)  focussed  the 
laser  beam  to  a  point  at  f]_,  introduced  a  knife  edge  at  the  light  source  image 
at  f£,  and  blocked  the  compensating  arm  of  the  interferometer.  The  quality  of 
the  resulting  pictures  (Fig.  cU,  C5  of  Ref.  l)  was  unsatisfactory  for  the 
present  study.  Igra  points  out  several  pertinent  factors:  l)  the  undiffused 
coherent  laser  light  and  the  very  small  source  size  caused  significant  dif¬ 
fraction  phenomena  to  appear  on  the  film  plate,  (Ref.  3*0;  2)  the  focussing 

of  the  powerful  laser  beam  caused  burning  at  the  knife  edge,  and  3)  attempting 
to  place  the  knife  edge  near  the  source  image  of  both  wavelengths  (f£  depends 
on  wavelength  and  |f^  -  f^o|  =  3mm)  gave  inferior  results  for  both.  All  these 
factors  contributed  to  the  detriment  of  the  schlieren  images. 

These  problems  were  avoided  in  the  present  study  by  focussing  with  a 
cylindrical  lens  the  powerful  laser  beam  onto  h  masked  glass  diffuser  to 
produce  a  sharply-defined  narrow  rectangular  light  source  of  less  coherent 
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monochromatic  light  at  t'^„  The  glass  diffuser  made  the  resistant  ambient  field 
of  the  teat-section  image  more  uniform  and  reproducible,  which  waB  of  great 
importance  in  Identifying  flow  phenomena  by  comparing  flow  pictures  with  no- 
tiow  pictures.  The  laser  beam  did  not  burn  the  glass 'diffuser  when  focussed 
to  the  narrowest  attainable  line,  but  the  large  focal  length  (175  mo)  of  the 
50  mm  x  60  mm  cylindrical  lens  may  have  prevented  precise  focusting.  The 
second  harmonic  of  the  laser  light  was  not  used  so  the  knife  edge  was  placed 
at  the  red  Xj-image  of  the  source  by  using  the  red  (X  *  6923  A)  continuous 
light  of  the  He-Ne  aiming  laser.  Retaining  or  removing  the  SHG  had  no  observable 
effect  on  the  schlieren  images. 

Since  the  visualized  phenomenon,  the  shock  front,  had  its  density 
gradient  in  the  x -direct ion,  that  is,  horizontal,  the  slit  source  was  oriented 
vertical,  parallel  to  the  knife  edge.  The  knife  edge  was  arranged  so  that  a 
positive  refractive  index  gradient  in  the  x-direction  would  cause  the  source 
image  at  f2  to  be  displaced  to  let  more  light  past  to  fall  upon  the  photo¬ 
graphic  plate.  The  camera  was  focussed  on  the  reference  wires  15  cm  from  the 
center  of  the  test  section  toward  the  camera  as  suggested  by  Hall  (Ref.  20). 

This  made  possible  accurate  spatial  interpretation  of  the  schlieren  images  by 
determining  rrom  the  real  wire  separations  and  their  sharp  images  on  each 
picture  the  magnify cati cn  of  that  picture.  Since  the  vertical  wires  were 
indeed  parallel  to  the  y-axis,  deviations  of  3hock-front  images  from  vertical 
could  be  easily  detected  and  analysed. 

Despite  the  above-mentioned  alsadvan* ages,  Tgra  used  no  diffuBer  in  tne 
laser  optics  in  order  to  generate  a  light  source  as  small  as  possible  to 
maximize  the  sensitivity  of  the  schlieren  system.  In  the  present  study  it  was 
felt  that  any  sensrtibity  1  36t  by  using  a  diffused  slit-source  instead  of  sin 
undiffused  point-source  would  be  gained  in  the  increased  quality  and 
uniformity  of  the  resulting  picturee.  A  maximum  BenBitivity  was  necessary  to 
obtain  good  schlieren  pictures  because  of  the  quite  low  densities  involved  in 
producing  strong  shocks  in  shock  tubes,  arid  because  of  the  low  energy  of  the 
normal -lasing  pulses  used. 

In  a  conventional  schlieren  system  the  knife  edge  is  adjusted  to  let 
about  half  the  light  from  the  rectangular  light  source  pass  by  to  fall  upon  the 
photographic  plate,  thereby  creating  an  ambient  or  "no-flow”  field.  Density 
gradients  of  one  sign,  which  are  normal  to  the  knife  edge,  cause  more  light  to 
pass  while  gradients  of  the  opposite  sign  cause  less  light  to  pass.  Thus  areas 
lighter  and  darker  than  the  ambient  field  are  seen  on  the  picture  and  are 
interpreted  as  corresponding  positive  and  negative  density  gradients  inside 
the  test  section. 

In  this  multiple-exposure  schlieren  Bystem,  however,  where  the  giant 
pulse  from  the  laser  is  much  stronger  than  the  subsequent  normal-lasing  pulses, 
an  ambient  field  of  significant  intensity  (about  half  total)  caused  by  the 
giant  pulse  would  overpower  the  images  of  the  translational  shock  front 
produced  by  the  weak  normal  pulses.  To  prevent  this,  the  knife  edge  was 
adjusted  to  block  off  more  of  the  source  image  so  that  the  ambient  field 
created  by  the  giant  pulse  was  very  weak  and  the  normal-lasing  pulse  Images 
would  show  up  clearly  through  it.  This  procedure  was  feasible  in  this  case 
because  the  flow  phenomenon  to  be  observed,  the  translational  shock  front,  had 
density  gradients  of  only  one  sign.  The  knife  edge  was  oriented  such  that 
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these  density  gradients  bent  the  light  past  the  knife  edge,  causing  the  shock 
front  to  appear  in  a  photograph  as  a  light  area  on  a  dark  ambient  field  (and 
opposite  on  the  film  plate).  Since  no  density  gradients  of  opposite  sign  were 
of  concern,  it  was  not  necessary  to  distinguish  areas  darker  than  the  ambient 
field  to  be  as  dark  as  was  thought  best.  Due  to  the  nature  of  the  light  pulses 
emitted  by  the  ruby  laser  operated  in  the  combined  giant-pulse/normal-lasing 
mode  (t  =  500  us),  this  modification  was  essential  to  the  success  of  the 
multiple -expo sure  technique. 

This  explains  why  the  schlieren  pictures  that  appear  in  this  experimental 
work  seem  to  be  unusually  dark.  They  generally  consist  of  a  background  that 
shows  the  elliptical  image  (owing  to  inclination  of  the  beam-splitters  and 
mirrors)  of  the  round  test-section  windows  and  the  reference  wires,  and  several 
brighter  shock-front  images.  The  number  and  brightness  of  the  images  depended 
upon  the  unpredictable  energies  contained  in  the  normal-lasing  pulses  described 
in  Sec.  2.2.1.  The  faint  horizontal  streaking  found  in  many  of  the  darker 
schlieren  pictures  was  caused  by  the  grain  in  the  glass  diffuser  when  the  source 
width  was  very  small,  but  this  nonuniformity  in  the  background  was  identical  in 
flow  and  no-flow  pictures,  and  was  not  an  obstacle  in  identifying  flow  disturban¬ 
ces.  Also,  the  shock  front  images  which  were  thin.,  vertical  lines  were 
perpendicular  to  these  streaks  and  the  former  were  easily  perceived  despite  the 
presence  of  the  latter.  A  very-fine-grain  diffuser  would  improve  future 
schlieren  pictures  obtained  in  this  manner. 

The  reason  for  maximizing  the  sensitivity  of  the  schlieren  system  and  for 
maintaining  a  very  faint  background  field  was  to  obtain  a  usable  contrast 
between  the  giant-pulse  ambient  field  and  the  normal-lasing  images  of  the  shock 
front  in  order  to  successfully  observe  the  shock-front  geometry  of  an  ionizing 
shock  wave  at  different  times.  This  contrast  was  enhanced  in  both  the  photo¬ 
graphic  and  printing  process.  Kodak,  high-density  B  &  W  negative  plates  were 
exposed  to  the  schlieren  image.  Prints  from  these  negatives  were  done  on 
Kodak  F-5  high-contrast  paper  and  given  a  high-gloss  finish.  The  plates  were 
capable  of  representing  a  very  wide  range  of  exposure.  On  the  other  hand,  the 
high-contrast  paper  could  represent  only  a  fraction  of  the  plate's  total  range, 
expanding  this  limited  range  to  the  complete  set  of  tones  of  the  paper.  The 
middle  of  this  fractional  range  was  determined  by  the  exposure  to  which  the 
paper  was  subjected  in  printing.  In  practice,  this  meant  that  in  spite  of  the 
wide  variation  of  absolute  exposure  on  the  film  plate  of  the  shock  images  and 
background  due  to  changes  in  both  controllable  and  random  optical  parameters 
(see  Appendix  III),  the  contrast  between  them  could  be  successfully  optimized 
by  varying  the  exposure  given  to  the  printing  paper.  Flow  and  no-flow  plates 
and  prints  were  processed  identically  to  minimize  differences  not  associated  with 
flow  phenomena,  although  differences  in  plate  exposure  did  occur  because  of 
slight  random  variations  in  laser  output  energy. 

Simultaneous  monochromatic  dual -wavelength  schlieren  pictures  of  the 
instabilities  studied  here  would  have  provided  information  on  the  refractive 
index  gradients  in  the  test  section  at  wavelengths  of  X^  and  Xg.  Since  the 
magnitude  of  the  electron  refractive  index  is  greater  at  X^  than  at  X2, 
information  about  the  presence  of  mass-density  gradients  and  electron-number- 
density  gradients  separately  in  the  disturbed  areas  behind  the  curved  translational 
shock  front  could  be  obtained.  This  is  practical  because  of  the  appreciable 
refractive  index  associated  with  the  electron-number  density  (five  times  as 
great  in  magnitude  than  that  of  the  mass  density  at  wavelength  X^  and  the 
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large  electron-number  densities  (up  to  2.5  x  10^-T  cm- 3)  associated  with  this 
phenomenon.  The  information  might  have  led  to  a  better  understanding  of  the 
possible  ion-acoustic  behaviour  of  the  plasma  behind  the  shock  as  related  to 
the  observed  geometry.  However,  there  are  several  disadvantages  to  this  method. 
The  presence  of  the  SHG  attenuates  the  laser  output  at  wavelength  X^,  resulting 
in  light  of  wavelength  X2  of  only  about  10/K  of  the  original  energy.  The  beam 
splitter  in  the  camera  attenuates  both  beams  again  by  at  least  one  half.  The 
fact  that  the  source  image  position  f£  is  not  the  same  for  both  and  X2  means 
that  the  knife  edge  must  be  placed  to  give  a  compromise  in  picture  quality.  All 
these  detract  from  the  sensitivity  of  the  system.  These  attenuations  would  most 
likely  prevent  the  appearance  of  the  fainter  images  produced  by  the  weaker 
normal-lasing  light  pulses,  certainly  for  the  second  harmonic,  if  the  multiple- 
exposure  mode  (tp  =  500  ys )  were  used.  However,  to  obtain  one  acceptable  image 
for  each  wavelength,  the  giant-pulse  mode  (tp  =  900  ps)  must  be  used,  at  the 
expense  of  any  information  about  the  time  behaviour  of  the  shock  front 
geometry. 

Since  it  was  this  time  behaviour  information  and  hence  the  multiple- 
exposure  method  that  was  desired,  the  second  harmonic  Xg  was  not  used.  The  SHG 
and  the  camera  beam  splitter  was  removed  to  prevent  any  attenuation  associated 
with  them.  The  knife  edge  was  placed  precisely  at  the  source  image  at  f2  since 
only  X-^  was  present.  All  these  measures  served  to  increase  the  source-image 
intensity  at  the  knife  edge  and  hence  increased  the  sensitivity  of  the  schlieren 
system  in  order  to  observe  the  images  produced  by  the  low-energy  normal-lasing 
pulses  when  tp  -  500  ps.  Thus,  knowledge  of  the  time  behaviour  of  the  shock 
front  was  gained  at  the  expense  of  observing  any  separable  density  and 
electron-number-density  gradients  in  the  complex  region  behind  the  non-planar 
shock  front, 

2.U  Shock-Wave  Signature  Technique 

The  technique  of  visualizing  the  pressure  pulses  and  slipstreams  associated 
with  detonation  waves  by  observing  the  imprints  left  on  a  film-coated  wall 
(Ref.  27)  was  attempted  here.  It  was  hoped  that  in  this  way  any  local  pressure 
pulses  or  slipstreams  caused  by  shock  intersections  linked  with  the  non-planar 
shock- front  geometry  would  appear,  that  is,  leave  a  "signature",  if  they  moved 
transversely  across  the  sidewall  of  the  shock  tube. 

Each  of  the  two  removable  test-section  window  frames  is  a  rectangle  18  cm 
(7")  high  and  38  cm  (15")  wide,  containing  in  the  center  the  window  which  is 
20  cm  (8")  in  diameter,  and  leaving  an  open  area  on  each  side  18  cm  high  and 
9  cm  (3.5")  wide.  The  downstream  area  of  each  was  chosen  for  the  coating  to 
avoid  any  damage  to  the  interferometric-quality  windows.  Due  to  the  high  gas 
speeds  and  very  high  temperatures  behind  the  shock  wave  (over  10,000°K)  it  was 
felt  that  no  Mylar  or  tape  should  be  left  in  the  tube  [Oppenheim  (Ref.  31) 
coated  sheets  of  thin  Mylar  which  were  taped  to  his  shock-tube  wall].  Thus,  a 
thin  film  of  soot  from  an  over-rich  acetylene  torch  was  applied  directly  to  these 
two  areas  as  evenly  as  possible.  Despite  the  utmost  care,  blacker  streaks 
made  by  the  thin  stream  emitted  by  the  torch  nozzle  were  evident.  The  area  on 
the  window  frame  which  was  removed  was  coated  most  successfully  and  photographed. 
The  corresponding  area  on  the  other  frame  which  was  not  removed  was  more  difficult 
to  reach  to  coat  and  photograph  because  of  the  presence  of  the  interferometer 
and  delicate  window  surface  only  centimeters  away. 
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The  window  frame  was  installed,  the  shock  wave  produced,  and  the  frame  was 
removed  and  photographed  again.  In  this  way  comparable  no-flow  and  flow  pictures 
of  the  soot-covered  wall  were  obtained. 


EXPERIMENTAL  RESULTS 


Performance  of  Multiple-! 


3 sure  Schlieren  System 


Twenty-three  runs  were  successfully  conducted  using  the  multiple-exposure 
schlieren  technique,  nineteen  of  which  resulted  in  pictures  showing  more  than 
one  image  of  the  translational  shock  front.  The  giant  pulse  from  the  Q- 
switched  ruby  laser  gave  a  good  schlieren  image  of  the  shock  each  time.  The 
clearness  or  absence  of  the  images  produced  by  the  subsequent  normal-lasing 
pulses  varied  from  run  to  run;  if  any  of  these  showed  up,  the  image 
corresponding  to  P2  was  certain  to  be  most  distinctive. 


To  examine  the  behaviour  and  effectiveness  of  this  technique  for  the 
purpose  of  multiple- expo sure  photography,  the  number  of  usable  shock-front 
images  in  each  picture  was  plotted  versus  the  pumping  time  tp  in  Fig.  17- 
from  data  in  Table  3.  A  shock- front  ii’age  was  considered  usable  if  its 
geometry  or  shape  could  be  observed  and  measured.  This  graph  shows  that  to 
obtain  the  multiple-exposure  objective  the  operational  parameter  tp  (defined 
in  Sec.  2.2.1)  is  to  be  kept  below  515-520  ys. 


The  brightness  of  the  normal-lasing  images  ranged  from  almost  the  same 
as  the  giant-pulse  image  to  barely  perceptible  but  still  usable.  Near 
tp  =  1+90  ys  the  image  produced  by  p2  was  good,  but  for  lower  t  all  the 
images  were  very  faint.  This  agrees  with  the  laser-pumping  described  in 
Sec.  2.2.1,  where  neither  the  energy  stored  in  the  ruby  crystal  was  due  to  popu- 
Irfiion  inversion  nor  the  power  available  for  normal  lasing  had  reached  its 


peak  before  the  maximum  flash-lamp  output  occurred  near  tp  =  U80  to  500  ys. 
On  the  other  hand,  after  about  tp  =  530  ys  the  flash-lamp  output  had  dropped 
but  the  energy  accumulated  in  the  crystal  had  increased,  and  only  a  bright 


giant-pulse  image  was  produced  upon  Q-switching.  This  was  indeed  observed. 


Thus  the  best  multiple-exposure  pictures  occurred,  as  expected,  for  tp 
between  1+70  and  520  ys.  All  runs  were  designed  for  t  =  500  ys,  but  the 
scatter  in  the  observed  values  is  due  to  1)  the  failure  to  predict  precisely 
the  diaphragm  breaking  pressure  and  hence  the  shock  speed,  2)  failure  to  fill 
the  channel  with  exactly  the  designed  initial  pressure  because  of  a  reading 
error  of  about  0.1  torr,  3)  variations  in  performance  of  the  shock  detectors, 
and  U)  uncertain  deceleration  of  the  shock.  Future  experiments  of  this  type 


should  design  for  about  tp  =  1+8 5  ys. 


Schlieren 


ses  of  Ionizing  Shock  Wave 


Runs  conducted  to  calibrate  the  schlieren  system  were  designed  for  lower 
shock  speeds  in  order  to  observe  the  behaviour  of  the  multiple-exposure  tech¬ 
nique  under  well-known  shock  conditions.  Figure  18  shows  a  multiple-exposure 
schlieren  picture  of  a  shock  wave  travelling  from  right  to  left  at  M_  *  2.99 
into  argon  at  an  initial  pressure  of  27.0  torr  and  temperature  of  299.1°K.  At 
this  Mach  number  the  only  (significant)  density  gradient  occurs  at  the  trans¬ 
lational  shock  front.  Each  shock-front  image  is  labelled  according  to  the  pulse 
from  the  laser,  operating  in  the  combined  giant-pulse /normal-lasing  mode 


VS 


(tp  =  527  Vs),  which  produced  it.  Thus  the  giant  pulse  created  hy  the 
Q- switch,  gave  the  image  labelled  1^.  The  first  normal-lasing  pulse  Pp 
produced  I2  some  time  later,  and  so  on.  Four  shock  images  appear  clearly  in 
Fig.  18,  and  much  fainter  ones  can  be  found  on  the  original  print.  This  is 
consistent  with  the  typical  laser  output  trace  in  Fig.  15 .  As  with  all  others 
of  its  type,  this  picture  was  printed  on  high-contrast  photographic  paper  so 
as  to  maximize  the  contrast  between  the  shock  images  and  the  background; 
otherwise  a  fainter  shock  image  such  as  Ig  is  indistinguishable  from  the 
background. 

The  shock  waves  appear  as  a  straight  white  line  on  a  darker  background 
because  of  the  orientation  of  the  knife  edge  as  discussed  in  Sec.  2.3*3.  The 
thickness  of  the  schlieren  image  of  the  shock  front  reflects  not  the  actual 
thickness  but  instead  optical  parameters,  such  as  the  source  size  (width)  and 
the  degree  to  which  the  shoek  surface  is  parallel  to  the  z-axis.  Note  the 
darkness  of  the  background  or  ambient  field  illuminated  chiefly  by  P-^,  and 
how  little  actual  contrast  there  is  between  Ig  and  this  field  which,  if  any 
lighter,  would  have  obscured  it.  This  shows  that  the  practice  of  cutting  out 
most  of  the  source  light  at  the  knife  edge  was  Justified  for  this  technique. 

Most  of  the  normal-lasing  images  were  in  fact  near  this  broaderline  contrast. 

In  addition  to  the  schlieren  image,  a  slight  shadow  graph  effect  appears 
at  each  shock  image  as  a  vea«y  thin  dark  line  to  the  left  of  and  immediately 
adjacent  to  the  thicker  white  line,  a  higher  density  behind  the  travelling  shock 
wave  as  expected.  A  v  inverse  shadowgraph  effect  is  caused  by  the  vertical  wire 
wv  but  not  by  the  horizontal  wire  w^  whose  shore  dimension  is  parallel  to  the 
knife  edge.  Compare  the  straightness  of  Ij_  to  Ij,  with  wv  which  is  nearly 
perfectly  vertical  and  serves  in  defining  the  y-direction.  The  wire  wg  is 
parallel  to  the  shock-tube  bottom  wall  which  is  visible  in  each  picture,  and 
serves  in  defining  the  x-direction. 

The  shock  velocity  in  this  case  is  Us  =  963  m/s  and  the  distances  between 
the  shock  images  can  be  measured  accurately.  Thus  the  times  between  the 
emission  of  the  giant  pulse  p;  and  the  three  pulses  P2,  P3,  and  P^  to  create 
the  four  images  of  the  travelling  shock  wave,  as  shown,  were  7-36  ys,  30.8ys, 
and  58.7  ys,  respectively.  It  was  soon  discovered  that  the  time  of  7.36  ys 
between  the  emission  of  the  first  two  light  pulses  was  very  typical,  giving 
only  a  horizontal  displacement  in  shock  images  at  Mach  3  of  7.1  mm;  but  at 
Mach  15,  where  ionization  occurs,  one  of  3.5^  cm. 

Figure  19  is  a  multiple-exposure  schlieren  picture  of  a  shock  wave 
travelling  from  right  to  left  at  Mach  6.38  into  argon  at  a  pressure  of  7.^5 
torr  and  a  temperature  of  299* 0°K.  The  shock-wave  images  1^  -  1^  are 
produced  by  the  successive  light  pulses  -  P^  from  the  laser  operated  with 
tp  =  U89  ys.  The  shock  velocity  is  2.02  x  103  m/s,  giving  the  time  between  the 
emissions  of  P-^  and  P2  as  6.9**  Vs,  again  near  7  Vs.  The  times  between  the  images 
other  than  I1  and  I2  are  not  reproducible  amongst  experiments,  as  can  be  seen 
by  comparing  this  picture  with  Fig.  18.  The  post-shock  temperature  T-  is 
3950°K  and  very  little  ionization  occurs  (x  <  10“°) .  Hence,  the  shock  wave  is 
observed  in  the  picture  to  be  perfectly  straight.  For. Mach  3  and  6.2,  removing  the 
Xj  band-pass  filter  which  screens  out  broadband  radiation  has  little  effect  since 
there  is  no  detectable  radiation  fYom  the  hot  gas  at  this  temperature.  In  these 
cases  the  laser/spark  light  source  described  in  Sec.  2.2.2  could  be  used,  but,  as 
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shown  in  Figs.  18  and  19,  no  non-planar  shock  geometries  appear. 

Figure  20a  is  a  multiple-exposure  schlieren  record  of  an  ionizing  shock  wave 
moving  from  right  to  left  at  Ms  »  11.5  into  quiescent  argon  at  a  pressure  of 
9.98  torr  and  a  temperature  of  299-2°K  (Run  #2  in  Tables  2  and  3).  There  are 
two  good  images  1^  and  I2  caused  by  and  Pg,  and  two  faint  ones  may  be  seen  at 
Ig  and  Ijj.  The  shock  images  are  thicker  here  than  in  Figs.  18  and  19.  This  is 
most  likely  due  to  the  shock  not  being  parallel  to  the  illuminating  light  beam, 
that  is,  the  z-direction  as  defined  in  Fig.  1.  The  greater  thicknesses  of  the 
succeeding  images  may  also  be  caused  by  normal-lasing  light  pulses  that  fail  to 
be  short  enough  to  "freeze"  the  shock  motion,  although  this  is  doubtful.  The 
width  of  the  rectangular  (slit)  light  source  which  is  located  at  the  interferometer 
input  focal  point  fi  was  reduced  significantly  from  that  for  Fig.  19.  The  increased 
schlieren  sensitivity  shows  the  vertical  wire  as  white  (schlieren  effect)  instead 
of  dark  (blocking  effect)  as  before. 

The  larger  image  thickness  does  not  prevent  us  from  detecting  and  describing 
the  geometry  of  each.  Using  the  vertical  wires  as  a  reference,  pairs  of 
vertical  lines  containing  the  shock  images  can  be  drawn  parallel  to  these  wires 
which  are  perpendicular  to  the  bottom  wall  evident  at  the  bottom  of  the  picture. 

(The  top  wall  does  not  appear  in  any  picture  since  both  cannot  be  viewed 
simultaneously  with  the  interferometer.)  Each  pair  of  lines  is  illustrated  in 
Fig.  20b  and  is  defined  as  follows:  The  left  one  1*  was  drawn  to  be  tangent 
to  the  shock  image  at  its  left-most  extent,  and  right  one  1'  was  drawn  to  be 
tangent  at  the  right-most  extent.  The  lines  in  Fig.  20b  for  1^  are  1J  and 
1^,  respectively.  Each  pair  of  lines  constructed  in  this  way  was  separated  by 
some  distance  in  the  x-direction  which  was  measured  and  divided  by  the  picture 
magnification  to  get  the  real  separation.  The  distance  between  these  two  tangents 
to  the  extreme  points  of  the  perturbed  shock  front  is  called  the  "total  perturbation" 
by  Lapworth  (Ref.  26).  Here,  we  define  the  amplitude  g  of  the  shock  front 
non-planarity  to  be  one  half  of  this  quantity.  Thus,  the  amplitudes  of  the  various 
images  of  the  same  shock  can  be  compared.  For  a  planar  shock,  of  course,  this 
amplitude  is  zero,  and  1*  and  1'  coincide.  The  shock  position  1  was  defined  as 
the  mean  position  of  1*  and  1'.  Thus  the  image  separations  are  defined,  and  the 
time  between  light  pulse  emissions  can  be  carried  out  as  before.  The  shock  speed 
for  Fig.  20  is  3.71  x  lO^m/s ,  so  the  time  delay  between  and  Pg  as  indicated  by 
1^  and  Ig  is  6.86  ys. 

The  shock  appears  to  be  perpendicular  to  the  bottom  wall  at  1J  and  approaching 
the  top  wall  in  the  same  manner  at  lj.  The  curve  traced  by  the  image  I,  is  very 
smooth,  with  no  kinks  (infinite  curvature)  that  indicate  the  presence  or  triple 
points  which  are  evident  in  pictures  of  multi-headed  gaseous  detonation  waves 
•(Refs.  32,33).  The  shape* of  the  curve  in  fact  suggests  a  sinusoidal  deviation 
from  a  straight  line,  consistent  with  the  observations  of  Bristow  (Ref,  2).  If 
the  origin  of  the  (x'-y')  plane  is  fixed  to  the  intersection  of  1-^  and  the 
bottom  wall,  then  the  shock  image  may  be  closely  described  by  the  locus  of  pointB 
satisfying 

x  =  g  cos(ky)  (3.1) 

1 

where  g  is  the  non-planarity  amplitude  described  above  and  k  *  tt/L  wheje  L  is 
the  height  of  the  shock  tube.  When  y  ranges  from  0  to  L,  the  angle  ky  goes 
through  it  radians  producing  one  half  of  a  complete  cosine  wave.  Since  any 
number  of  these  half-waves  can  be  connected  to  fill  the  length  L  and  still 
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satisfy  the  boundary  condition,  they  are  used  as  the  unit  to  describe  the  shock 
image  geometry.  The  number  of  half-waves  nh  is  the  number  of  these  geometrical 
units  detected  in  a  shock  image.  Thus  the  number  of  half-waves  in  image  1^  in 
Fig.  20  is  nh  =  1. 

A  curve  with  n.  of  these  units  between  y  =  0  and  y  =  L  is 


x  =  g  cos(nhky)  (3.2) 

+ 

The  parameter  nfa  is  also  equal  to  the  number  of  times  this  curve  crosses  x  =  0. 

If  an  actual  shock  image  I  at  1  is  close  to  sinusoidal,  then  counting  the 
intersections  of  1  with  the  wavy  image  I  to  find  nh  is  easier  in  practise  than 
other  methods  involving  points  of  tangency  mentioned  in  the  definitions  of  1* 
and  1'.  Identifying  the  exact  geometry  of  each  image  of  the  same  shock  and 
subsequently  determining  the  change  over  time  is  also  done  more  reliably  by 
examining  the  I  -  1  intersections.  Uncertainties  in  the  value  of  the  interger 
n,  are  due  to  the  non-sinusoidal  shape  of  I  where  intersections  with  the 
straight  line  1  are  difficult  to  determine.  The  maximum  error  in  this  study 
is  ±1,  but  most  images  have  a  well-defined  value  for  n^. 

Figure  21a  is  the  multiple-exposure  schlieren  picture  of  an  ionizing  shock 
wave  travelling  from  right  to  left  at  Ms  =  l4 . 1*  into  argon  at  a  pressure  of 
Uill  torr  and  a  temperature  of  297. 0°K  (Run  #8).  A  value  of  50U  ys  was 
observed  for  t  .  The  giant  pulse  P-j_  gave  a  bright  image  I1  (compare  with  1^ 
in  the  three  previous  pictures),  but  no  other  images  are  detectable.  Since  the 
shock  speed  was  U.62  x  lO^m/s  the  consistent  time  delay  of  about  7  US  between 
and  P2  would  cause  I 2  to  appear  near  the  place  marked  ”12".  The  absence  of 
any  secondary  images  is  probably  due  to  the  weakness  of  the  normal-lasing 
pulses  for  tp  over  500  ys.  The  value  of  nh  =  3,  and  g  is  measured  as  1.6  mm. 

Behind  the  shock  image  there  is  a  system  of  intersecting  lines  that 
seem  to  indicate  transverse-moving  and  colliding  shocks.  Disturbances 
behind  a  steady  shock  in  the  shape  of  a  half-wave  would  be  caused  by  rarefaction 
and  shock  (compression)  waves.  These  components  are  shown  in  the  drawing 
Fig.  21b  for  clarity. 

A  vertical  band  D-,  about  1.5  cm  wide,  that  is  darker  than  the  background, 
appears  about  U  cm  behind  the  shock  image,  where  the  lines  W-^  are  found.  This 
feature  is  conspicuously  absent  in  the  no-flow  picture.  This  dark  area  is  most 
likely  due  to  the  large  electron  number  density  gradient  dne/dx  (electron 
cascade)  Just  before  ne  reaches  its  maximum  value  ne,E  at  the  point  of  quasi¬ 
equilibrium  behind  the  strong  shock.  Since  the  electrons  make  a  negative 
contribution  to  the  refractivity  of  the  gas,  a  positive  dne/dx  would  cause  a 
negative  refractive  index  gradient,  leaving  a  dark  area  on  the  schlieren  record. 
This  explanation  coincides  with  the  fringe  shifts  observed  in  the  monochromatic 
interferogram  in  Fig.  22  of  a  similar  ionizing  shock  in  argon  where  Mg  *  lh.77, 

Pl  =  3.95  torr,  and  T^  *  297»6°K.  The  large  dne/dx  is  indicated  by  a  large 
sudden  fringe  displacement  in  a  direction  opposite  to  that  occurring  at  the 
translational  shock  front.  Note  that  the  lines  W-j_  are  much  brighter  inside  the 
band  Dj.  This  signifies  that  these  transverse  waves  or  shocks  are  intensified 
in  this  region.  This  effect  can  also  be  seen  in  the  interferogram  of  Fig.  22, 
in  which  the  fringes  are  disturbed  most  in  this  region  of  sudden  fringe  shift, 
caused  by  the  electron-cascade  front  where  n@  achieves  its  equilibrium  value. 
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Figure  23  is  a  multiple-exposure  schlieren  record  of  an  ionizing  shock  wave 
moving  from  right  to  left  at  M_  =  l4.5  into  argon  at  a  pressure  of  3.18  torr  and  a 
temperature  of  299. 0°K  (Run  #7J.  The  most  notable  features  of  this  picture  are 
the  large  amplitude  of  I^,  the  closeness  in  shape  of  1^  to  an  ideal  cosine  curve 
with  njj  =  2,  and  the  uniform  image  width  of  I3.  The  amplitude  g  of  1^  is  3.2  mm. 
and  the  image  is  sharp  because  of  the  very  short  duration  of  the  giant  pulse  P^. 
Image  Ig  is  faint  but  of  comparable  width,  with  the  same  n^  but  different  g. 

Image  I3  appears  to  have  been  exposed  by  a  light  pulse  of  about  0.7  ys  in 
duration.  Despite  this,  a  value  for  g3  and  nh  can  be  determined  without  difficulty. 
Since  the  shock  speed  was  4.66  x  lO^m/s,  the  separation  of  I^  and  Ig  gives  a 
delay  time  between  and  P2  of  6.57  ys . 

Figure  24  is  a  multiple-exposure  schlieren  picture  of  an  ionizing  shock 
wave  travelling  from  right  to  left  at  Ms  *  12.3  into  argon  at  a  pressure  of 
17. 04  torr  and  a  temperature  of  298. 5°K  (Run  #13).  Since  the  shock  speed  was 
396  x  lO^m/s,  the  time  between  P^  and  P2  was  6.64  ysec.  In  this  case  all  three 
shock  images  have  a  large  thickness,  even  I-,  which  was  exposed  by  the  short- 
duration  P-p  Thus  this  image  thickness  must  be  done  to  a  three-dimensional  effect 
where  the  shock  viewed  exposes  its  nonplanarity  in  two  directions  as  well. 

Other  clues  to  this  phenomenon  are  the  decreased  image  thickness  of  1^  near  the 
top  of  the  picture  and  the  increased  image  brightness  at  the  edges  of  the  image 
where  the  gradients  are  greatest.  This  is  exactly  what  is  seen  in  the  inter- 
ferogram  in  Fig.  25.  Determining  a  value  for  nh  for  each  shock  image  is 
possible,  but  identifying  g  is  somewhat  ambiguous.  Equation  3.2  describes  the 
locus  of  points  of  the  shock  translational  front  in  the  shock  tube  only  when  the 
shock  is  two-dimensional,  that  is,  independent  of  z,  as  in  1^  of  Fig.  23.  The 
non-planarity  in  the  z-direction  is  well  shown  by  Ij  in  Fig.  25.  Here  Ms  *  12.7, 
p-^  =  11.52  torr,  and  T-^  =  298. 1°K  (Run  #10);  t  =  511  ys  and  only  one  image  I2 
appears.  The  edges  of  1^  are  very  well-defined,  and  seem  to  constitute  two 
separate  shocks  with  nh  =  2  and  different  amplitudes.  The  wave  system  behind 
the  curve  shock  front  is  similar  to  that  observed  in  Fig.  22. 

Figure  26  is  a  multiple-exposure  schlieren  record  of  an  ionizing  shock  wave 
moving  from  right  to  left  at  Ms  =  16.I  into  quiescent  argon  at  5.79  torr  and 
299. 4°K  (Run  #20).  Two  usable  images  I^  and  I2  of  the  translational  shock  front 
appear,  each  with  n^  =  6.”  Each  has  behind  it  a  complex  system  of  disturbances 
(W^  and  W2)  intimatedly  linked  with  it.  The  shock  has  moved  forward  far 
enough  that  there  is  insignificant  overlap  between  the  schlieren  images  of  the 
complex  reacting  shock  produced  by  light  pulses  P-j_  and  P2  which  were  separated 
by  7.37  ys.  The  shape  of  I.  is  descri-ed  by  Eq.  3.2  quite  well  with  n^  =  6,  and 
the  disturbances  behind,  which  appear  to  be  converging  compression  waves  and 
intersecting  transverse  shocks,  are  related  to  I^  in  the  same  manner  as  in 
Fig.:'22  and  24.  Again,  the  basic  geometric  unit  seems  to  be  the  half-wave, 
and  the  compression-shock  waves  are  observed  for  each  half-wave  in  I^.  Each 
shock  in  turn  intersects  with  its  neighbour  to  produce  the  characteristic 
patterns  so  apparent  here  in  both  and  The  disturbances  W2  are  related 

in  exactly  the  same  way  to  the  geometry  of  Ig.  The  images  of  these  disturb¬ 
ances  are  remarkedly  well-defined  and  Bharp.  The  cross-hatching  observed  here 
bears  a  resemblance  to  the  pattern  left  by  a  multi-headed  gas-detonation  wave 
on  a  soot-covered  wall  (Refs.  31,32). 

Figure  27  is  a  multiple-exposure  schlieren  picture  of  an  ionizing  shock 
wave  travelling  from  right  to  left  at  Ms  =  14.5  into  quiescent  argon  at  11.35 
torr  and  300.1°K(Run  #21).  Five  distinct  shock“front  images  are  visible,  and 


some  parts  of  the  disturbances  behind  each  can  be  seen  as  well.  I,  is  clear  and 
shovs  some  three-dimensional  effects.  The  intensification  of  the  transverse 
shocks  in  the  electron  cascade  region  is  clearly  seen.  This  may  be  coupled  to 
the  transition  from  a  compression  wave  to  a  shock  wave.  Images  I2  and  1 3  are 
very  close  together  and  are  not  discemibly  different,  appearing  as  a  "double” 
image.  This  cescription  applies  as  well  to  the  corresponding  disturbances  W2 
and  W3  behind.  Images  1^  and  1^  are  similar  to  the  previous  pair  but  are 
slightly  more  separated.  Only  one  part  of  each  of  and  is  observable. 

P-L  and  P2  were  separated  by  7-1*1*  ys. 

All  interferograms  and  schlieren  pictures  of  ionizing  shocks  in  argon  had 
-some  degree  of  curvature  in  the  shock  front.  It  should  be  noted  that  for  each 
one  the  shock  front  was  perpendicular  to  the  visible  bottom  wall  at  the  wall. 

Sometimes  this  was  achieved  in  only  the  last  few  millimeters  from  the  wall  as  a 
required  boundary  conditions.  Also,  no  kinks  (infinite  curvature)  or  evidence 
of  triple  points  was  found  in  the  schlieren  image  of  any  shock  front  observed  in 
this  study.  This  can  also  be  said  for  all  such  interferograms  obtained  at  this 
facility  except  for  a  few- where  the  fringe  spacing  makes  interpretation  of  the 
exact  position  of  the  translational  shock  front  difficult.  This  shows  that 
although  the  nonplanar  shock  fronts  observed  here  seem  similar  in  appearance 
to  the  multi-headed  gas-detonation  waves  described  and  illustrated  in  Figs.  5  and 
7  of  Ref.  32,  there  is  to  essential  difference  indicated  by  the  lack  of  triple 
points  in  our  cases. 

•3.3  Temporal  Behaviour'  ef  Mon-Planar  Shock  Front 

Single-exposure  interferograms  and  multi pie-expo Bure  schlieren  pictures  with 
only  one  visible  image  wouid  give  us  no  definite  information  about  the  time 
behaviour  of  nonplanar  shock  fronts.  Only  those  schlieren  pictures  with  more  than 
one  useful  image  are  considered  in  this  section.  The  sinusoidal  Bhapes  of  each 
shock  front  observed  in  all  the  interferograms  of  ionizing  shocks  in  argon 
obtained  previous  to  this  study  led  to  the  following  three  hypotheses  regarding 
their  time  behaviour: 

1)  The  shock  front  was  curved  but  quasi-stationary 
(stationary  in  the  shock- fixed  coordinates). 

2)  The  shock  front  behaved  like  a  standing  cosine  wave. 

3)  The  shock  front  behaved  like  a  travelling  cosine  wave. 

A  shock  wave  is  either  perpendicular  to  or  reflected  by  a  solid  wall  to  main¬ 
tain  an  overall  flow  deflection  of  zero  at  the  wall.  In  the  absence  of  any 
reflections,  hypotheses  (l)  tod  (2)  would  be  consistent  with  this  boundary  condition 
but,  strictly  speaking,  (3)  would  not.  No  triple  points,  Mach  stems,  or 
reflections  were  found  near  the  wall,  but  occasionally  a  shock  front  image  was 
oblique  for  all  but  a  few  millimeters  away  from  the  wall.  The  amplitude  of  a  standing 
wave  is  near  zero  for  a  fraction  of  the  time,  and  a  series  of  single-exposure 
random  pictures  should  show  some  in  which  the  amplitude  is  very  small.  However, 
this  was  not  the  case.  Consequently,  single-exposure  pictures  were  insufficient 
to  determine  the  propagation  mechanism,  as  realized  by  Bristow  (Ref.  l). 
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Equation  3.2  describes  the  idealized  shock- front  geometry  at  one  moment, 
in  shock- fixed  coordinates.  The  evolution  in  time  of  the  shock  front  for  the 
above  hypotheses  would  be  described  in  shock-fixed  coordinates  as 

1)  x"  =  g  cosfn^ky")  (3.3) 


2) 


3) 


for  a  quasi-stationary  sinusoidally  curved  shock  front; 
x'  =  g  cos(n^ky')coswgt 

for  a  standing  cosine  wave  where  T  =  — -  is  the 

s  U) 

s 

period  of  oscillation;  and 
x'  =  g  cosfn^ky'  +  oi^t) 

for  a  travelling  cosine  wave  whose  phase  velocity  is  V 

w 

whose  period  is  =  X/V^. 


(3.U) 


(3.5) 


and 


All  the  multiple-exposure  pictures  were  examined  with  these  three  models  in 
mind.  Observations  of  the  change  in  geometry  of  the  shock  front  with  time 
consisted  of  the  measurements  of  several  key  parameters,  suggested  by  the  above 
three  models,  of  the  two  shock-front  images  I,  and  I2  which  were  created  by  the 
ruby  laser  giant  pulse  and  the  first  normal-lasting  pulBe  P2  that  occurred 
about  7  Us  after  P,  .  These  measurements  are  listed  in  Table  3,  the  Run  numbers 
corresponding  to  those  in  Table  2.  Parameters  tp  and  h8  are  discussed  in  Sec.  3.1. 
The  time  At  between  laser  pulses  P^  and  P2  is  the  quotient  of  the  distance 
between  the  mean  positions  1^  and  I2  in  the  test  section  of  the  nonplanar  shock- 
front  images,  and  the  shock  velocity.  Except  for  Run  #10,  this  delay  time  is 
quite  consistent. 

The  relationship  between  pumping  time  tp  and  the  number  of  distinctive  shock 
images  ns  is  illustrated  in  Fig.  17.  These  findings  may  be  useful  to  future 
researchers  considering  the  ruby  laser  used  in  this  combined  mode  as  a  multiple- 
exposure  light  source.  The  parameter  g  is  described  in  Eq.  3.1,  and  the  larger 
of  the  two  for  1^  and  Ip  is  given.  The  number  of  half-waves  n.  determines  the 
wavelength  X  of  each.  Note  that  each  shock  image  for  a  particular  run  had  the 
same  value  of  njj.  Hence,  the  value  of  n^  for  any  run  is  well-defined  and  is 
applicable  to  all  shock  images  in  the  photograph  obtained  from  that  run. 

The  pairs  of  shock  images  could  be  divided  roughly  into  three  groups:  those 
exhibiting  travelling-wave  behaviour,  those  exhibiting  standing -wave  behaviour, 
and  those  that  show  neither.  In  the  first  group,  the  images  1^  and  I2  are 
essentially  similar  but  displaced  vertically,  suggesting  a  change  in  phase  of  a 
travelling  wave.  The  measured  vertical  displacement  between  images  divided  by 
At  is  the  phase  velocity.  The  period  T+  in  Table  3  in  this  case  is  the  time  for 
the  shock  front  to  be  displaced  vertically  by  one  wavelength  from  1^;  the  wave 
type  is  denoted  as  T.  An  example  of  this  type  appears  in  Fig.  27,  which 
corresponds  to  Run  #21.  The  two  images  appear  to  be  quite  smooth  and  sinusoidal 
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in  shape,  but  appear  to  be  reversed  in  phase.  However,  examination  of  1^  reveals 
that  the  antinode  at  C  is  displaced  from  the  antinode  E  on  I»,  and  hence  must  have 
come  from  the  antinode  A  on  1^.  A  comparison  of  nodes  and  other  features  of  1^ 
with  those  of  In  and  Ijj  gives  a  similar  result.  The  vertical  displacement  Ay' 
between  and  I2  is  3.U7  cm  ±135>,  the  phase  velocity  is 

Vt  =  Ay'/At  =  3.U7/7.UU  =  0.1*7.  cm/  Us 
and  the  period  is 

Tt  =  X/Vt  =  7. 1/0. 1*7  =  15.1  Us 

In  four  runs  among  a  total  of  23  runs  performed  the  resulting  pair  of  shock- 
front  schlieren  images  suggested  a  travelling  sinusoidal  wave  as  opposed  to  a 
standing  wave,  that  is,  a  definite  phase  displacement  as  well  as  equal  amplitude 
was  observed.  It  should  be  noted  that  the  system  of  intersecting  shocks  behind 
the  curved  shock  front  (e.g.,  in  Fig.  27)  was  consistent  with  the  latter's 
geometry  for  each  shock  image.  For  each  I2  displaced  in  phase  from  Ii ,  the  system 
W2  was  moved  in  the  same  way  from  W]_.  This  is  further  proof  that  thecurved  shock 
front  and  the  regular  pattern  of  disturbances  behind  it  are  intimately  linked.  In 
each  case  observed,  the  phase  displacement  between  1^  and  I2  was  downward.  This  is 
reflected  in  the  plus  sign  which  appears  in  Eq.  3*5. 

In  the  second  group  the' two  shock  images  I^  and  I2  have  different  amplitudes 
but  seem  to  have  the  same  phase,  i.e.,  a  node  or  antinode  would  occur  at  the  same 
height  (y  value),  su-gesting  the  evolution  of  a  standing  wave.  Therefore, 
information  about  the  change  with  time  of  this  model  comes  from  the  two  distinct 
values  of  g  only.  In  all  cases  there  was  a  larger  value  and  a  smaller  value.  In 
order  to  obtain  a  rough  idea  of  the  period  in  this  mode  of  motion,  the  larger 
value  is  assigned  to  the  maximum  amplitude,  so  that  (x1 )mwv  =  g  in  Eq.  3.1*.  Since 
the  time  between  images  is  At,  the  angular  velocity  is  ios  =  A0/At,  and  the  period 

quoted  is  T  =  —  =  Atf£. 

s  0)  A0 
s 

In  Fig.  20  (Run  §2  in  Tables  2  and  3)  the  images  1^  and  I2  have  g  measurements 
of  1.8  mm  and  1.1  mm.  If  this  shock  was  a  standing  wave  with  amplitude  1.8  mm, 
then  the  angular  displacement  between  them  is  AS  =  cos~l 

o-°-» 

radians.  Since  this  change  occurred  in  t  =  6.86  ys,  the  period  is  Ts  =  6.86  ys 
x(i*/A0)  =  U8  ys.  The  fourth  image  1^,  which  has  g  =  0,  can  be  used  to  chick 
this  result.  It  is  U.67  cm  from  1^,  and  consequently  was  exposed  12.6  ys  after 
I1>  This  change  in  angle  is  one  quarter  of  a  period,  so  one  period  is  U  x  12.6  ys 
=  50.  U  ys,  a  time  quite  consistent  with  the  preceding  result.  In 
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this  way  the  presence  of  additional  images  again  increases  the  certainty  of 
the  applicability  of  the  model. 


There  are  some  multiple-exposure  schlieren  pictures  which  definitly 
show  a  travelling  wave  and  some  in  which  a  standing  wave  is  indicated. 

Others  do  not  show  definite  evidence  of  either  type  and  are  marked  in 
Table  3  by  no  wave-type  designation.  These  pictures  either  show  no 
discernible  difference  between  1^_  and  I2,  or  they  show  a  difference  accountable 
by  both  models.  Doubtful  decisions  are  accompanied  by  a  question  mark. 

3.1*  Shock-Wave  Signature  Records 

As  described  in  Sec.  2.U,  a  method  for  recording  the  possible  trajectory 
of  a  shock-front  triple  point  on  a  layer  of  black  soot  on  the  tube  wall  was 
implemented  here.  Great  care  was  taken  to  obtain  a  uniform  soot  layer,  but 
broad  black  streaks  are  seen  in  pre-run  photos  of  the  blackened  area.’  These 
photos  were  compared  with  post-run  pictures  of  the  area  to  detect  systematic 
transverse  traces.  It  was  assumed  that  such  traces  would  be  consistent  with 
the  geometry  of  the  same  shock  wave  observed  by  normal  photographic  means  in 
the  window  located  immediately  upstream. 

A  high  driver  pressure  p^  was  necessary  to  produce  conditions  suitable  for 
shock- front  oscillations.  This  required  the  combustion  of  a  large  amount  of 
O2  and  H2»  which  resulted  in  a  great  quantity  of  water  vapour  in  the  shock  tube 
immediately  after  an  experiment.  Condensation  caused  the  interior  walls  to  be 
coated  with  water.  The  imprints  were  preserved  by  opening  the  tube  and  allowing 
the  inside  to  dry  slowly. 

Three  experiments  were  conducted  using  this  technique.  An  inspection  of 
the  blackened  areas  after  each  run  revealed  no  markings  that  could  be  interpreted 
as  a  triple-point  trajectory  trace.  The  chief  character  of  these  areas  was 
that  of  a  multitude  of  horizontal  steaks  supposedly  created  by  the  passage  of 
the  wet  expanding  driver  gas. 

It  is  concluded  that  this  visualization  technique  was  not  appropriate  for 
oscillating  ionizing  shocks  generated  in  this  manner.  Another  reason  for  the 
lack  of  triple-point  traces  may  simply  be  due  to  the  absence  of  any  triple 
points.  This  conjecture  agrees  with  the  observation  that  no  triple  points  were 
detected  in  interfdrograms  of  the  shock  front  taken  in  the  same  three 
experiments.  Also,  no  triple  points  were  observed  in  any  schlieren  or  interfero¬ 
metric  records  taken  of  this  oscillating  ionizing  shock  front. 

h.  Discussion  and  Conclusions 


The  problem  studied  here  could  be  classified  under  several  categories.  The 
shape  and  size  of  the  nonplanarity  of  the  translational  shock  front  suggest  that 
a  stability  analysis  might  be  in  order.  The  hydrodynamic  effects  of  a 
sinusoidally  curved  shock  front  on  the  relaxation  region  behind  could  be 
considered.  One  could  also  apply  the  principles  of  plasma  acoustics  to  the  two- 
temperature  ionizing  gas  in  order  to  predict  the  origin  and  type  of  oscillations. 
Each  of  these  approaches  will  be  considered,  relating  the  phenomenon  observed  to 
the  most  relevant  literature. 


Many  theorists  have  concerned  themselves  with  the  stability  of  a  plane 
shock  wave  faef.  35-*3).  Of  the  several  different  types  of  instability 
considered,  the  one  applicable  to  this  phenomenon  was  called  by  Freeman 
(Ref.  35)  "stability  jf  form".  He  looked  at  the  shock  wave  resulting  from 
moving  a  corrugated  piston  impulsively  into  a  fluid.  Small  deviations  of  the 
piston  surface  from  planar  allowed  a  linearized  theory  which  showed  that  the 
shock's  corresponding  perturbations  from  a  planar  form  decayed  in  time  like 
t~T  in  the  strong  shock  limit,  where  t  is  the  time  between  the  formation  of  the 
shock  and  the  measurement  of  the  perturbations.  This  technique  was  applied  to 
the  mere  experimentally  tenable  geometry  of  a  weak  shock  passing  over  thin 
wedges  on  the  walls  of  a  shock  tube,  with  similar  results  (Ref.  36).  The 
corresponding  experiments  were  done  by  Lapworth  (Ref.  26)  whose  composite 
picture  (see  Sec.  13)  depicts  a  shock  which  seems  to  undergo  several 
oscillations  somewhat  like  a  standing  wave,  with  two  half-wavelengths  (n^  =  2 
in  Eq.  3.4)  contained  in  the  5-875  in  high  channel.  In  fact,  a  model  similar 
to  Eq.  3.U  with  g  decreasing  with  time  was  used  to  describe  the  observed  shock. 
Closer  inspection  yields  some  important  differences,  however,  between  this  and 
the  presently  studied  phenomenon.  Lapworth' s  shocks  have  pronounced  kinks  and 
reflected  (or  transverse)  shocks  projecting  backwards  from  these  kinks.  No 
kinks  or  triple  points  were  observed  in  the  shocks  of  this  study.  Lapworth' s 
shock  waves  were  in  the  perfect  gas  regime  for  air,  with  Mach  numbers  less  than 
2,  as  opposed  to  the  strong  ionizing  shocks  in  argon  used  here.  Also,  he 
reported  no  travelling -wave  behaviour.  With  these  disparities  in  mind,  it  seems 
unlikely  that  the  two  instabilities  utilize  the  same  propagation  mechanisms  or 
sources,  despite  their  superficially  similar  appearance  and  time  behaviour. 

Thus,  a  stability-cf-form  analysis  as  applied  to  the  present  problem  would  be  of 
limited  usefulness- 


Other  researchers  considered  the  passage  of  a  shock  wave  through  a  medium 
with  an  arbitrary  equation  of  state,  examining  the  kinds  of  state  equations  and 
small  disturbances  that  v;.uid  lead  to  instability,  A  fundamental  paper  by 
D' Iakov  was  revised  and  clarified  in  an  illuminating  paper  by  Swan  and  Fowles 
(Ref.  38)  who  used  as  an  initial  perturbation  to  the  plane  shock  wave  a 
travelling  wave  equivalent  tz  that  expressed  in  Eq.  3.5.  (Remember,  Eq.  3.5 
is  in  a  sense  empirical,  since  it  derives  from  shocks  observed  in  this  study.) 
They  found  that  a  linear  analysis  resulted  in  stability  criteria  involving  the 
gradient of  the  Rankme-Hugoniot  curve,  the  Mach  number  M  of  the  fluid 
behind  the^shock  in  shock -fixed  coorinates,  and  the  slope  -j  of  the  Rayleigh 
line,  the  straight  line  on  the  p-v  plane  joining  the  pre-shock  state  (vp»Pr ) 
with  the  pest  shock  state  (v2,p2).  Here  v  is  the  specific  volume.  Fi&ur£  28 
illustrates  these  quantities,  including  the  mass  flux  j  through  the  shock. 

With  A  -  the  criterion  for  stability  of  the  shock  was  found  to  be 
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They  show  that,  if  A  violates  the  lower/upper  limit,  then  the  shock  could 
degenerate  into  two  waves  whose  directions  of  travel  are  similar /opposite. 


Griffiths  et  al  (Rtf.  39)  obtained  time-resolved  interferograms  of  very 
strong  shocks  in  argon  in  a  cylindrical  shock  tube.  A  shock  was  deemed  unstable 
if,  after  it.  reached  the  window,  the  shock  front  or  flow  behind  it  was  nonplanar 
but  distinct  from  the  contact  surface.  Unstable  shocks  were  observed  for 
shock  Mach  numbers  in  the  range  31  <  Ms  <  3*+,  which  coincided  with  velocities 
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such  that  the  first  ionization  of  argon  was  near  completion,  causing  (§^)h  to 
become  positive.  This  observed  stability  limit  was  related  to  inequality 

in  Eq.  U.l.  Their  data  suggested  the  different  upper  limit  stability  condition 
A  <  A' ,  where 

1  -  M^l  +  w./w  ) 

A’  =  - 3 - L-Z-  (U.2) 

1  -  M  (l  -  W1/w2) 

and  where  w^  and  W2  are  the  pre-  and  post-shock  gas  speeds  relative  to  the 
shock.  The  velocity  range  of  observed  instabilities  occurred  when  A  approached 
A'  closely.  The  equilibrium  post-shock  state  was  used  in  each  case  in  Ref.  39, 
but  no  comment  was  made  concerning  the  frozen  state  immediately  behind  the 
translational  shock  front,  or  the  two-temperature  plasma  in  the  ionization 
relaxation  region.  As  well,  no  attempt  was  made  to  link  the  non-plan ar  shock 
front  and  the  equilibrium  plasma  behind  it.  Unlike  the  present  study,  they 
observed  several  shocks  with  planar  translational  fronts  accompanied  by 
nonplanar  flow  behind. 

The  shock  waves  produced  in  this  study  were  not  near  completion  of  the 
first  degree  of  ionization.  That  is  X  <  0.2  (Tables  1  and  2).  The  maximum 
Mach  number  was  18.0,  far  short  of  the  observed  region  of  instability  of  Ref.  39. 
However,  their  shocks  with  disturbances  attributed  to  the  contact  surface,  as 
opposed  to  those  with  (^-)H  >  0,  were  placed  in  a  separate  class.  It  is 
interesting  to  note  that  Ref.  39  assigns  all  shocks  with  Mach  <  25  to  this 
class.  In  other  words,  no  stable  ionizing  shocks  in  argon  appeared  below 
Ms  =  25  to  their  minimum  of  Mg  =  18.  Whether  or  not  these  disturbances,  due 
to  the  contact  surface,  are  in  fact  similar  in  appearance  to  those  observed 
in  this  study,  cannot  be  seen.  Although  only  a  limited  range  of  Mach  numbers 
considered  by  the  two  investigations  were  in  common,  both  observed  that  ionizing 
shocks  in  argon  around  Ms  =  18  were  unstable.  A  more  definitive  look  at  this 
"contact  surface  disturbance"  would  be  in  order,  possibly  with  space-resolved 
interferograms .  It  is  safe  to  conclude  that  the  present  nonplanarities  cannot 
be  due  to  the  completion  of  the  first  degree  of  ionization  as  observed  in 
Ref.  39. 

It  should  be  noted  that  none  of  the  shock-wave  stability  theories, 
including  those  using  arbitrary  state  equations,  took  into  account  explicitly  a 
nonhomogeneous  or  a  nonequilibrium  post-shock  state,  such  as  the  two  temp¬ 
erature  ionization  relaxation  zone  behind  the  shock.  Whether  or  not  these 
effects  could  be  included  somehow  in  the  Hugoniot  curve  or  the  state  equation  is 
not  clear.  Nonequilibrium  studies  have  generally  been  in  the  areas  of  one¬ 
dimensional  shock  structure  (Ref.  2,3,11,23,25,1+1+)  and  acoustic  stability 
(Ref.  1+5, U6).  Also,  each  theory  assumed  shocks  and  perturbations  of  infinite 
extent,  so  that  sidewall  and  reflection  effects  were  ignored.  As  noted  before, 
when  a  travelling  wave  perturbation  was  observed  in  this  study  by  the  multiple- 
exposure  method,  reflections  of  the  perturbations  off  the  walls  were  not  evident, 
making  the  shock  front  conform  to  Eq.  3.5  more  closely  than  would  be  expected 
in  a  shock  tube  of  finite  dimension. 

The  hydrodynamic  consequences  of  a  sinusoidally  curved  shock  wave  can  be 
used  to  explain  the  compression  waves  immediately  behind  the  shock,  which 
converge  to  give  transverse  shocks.  These  intersect  with  their  neighbours  to 


give  the  system  of  waves  denoted  Wn  in  Figures  2  ,  25,  26,  and  27,  where,  for 
example,  corresponds  to  the  image  1-^  of  the  translational  shock  front.  In 
order  to  understand  this  wave  system,  consider  the  stationary  shock  wave  in  Fig. 
29,  where  the  amplitude  of  the  perturbation  has  been  exaggerated  for  clarity. 

The  gas  enters  from  the  left  and  passes  through  the  shock  at  an  angle  3,  where 
3  =  ^  for  a  normal  shock.  Since  tan  (3)  i  cot(S)  =  and 

3  =  cot-1[-gn^k  sin(n^ky) ]  (U.3) 

We  can  see  from  the  diagram  that  the  shock  is  normal  when  cos(n^ky)  =  ±1, 
coinciding  with  sinCn^kv)  =  0,  and  giving  3  -  cot~^(0)  =  5-,  as  expected.  The 
quantity  gk  is  dimensionless,  representing  the  relative  magnitude  of  the 
perturbation.  The  maximum  deviation  of  the  shock  front  from  normal  is  found 
for  the  one  in  Fig.  26,  for  example,  where  n^  =  6  and  g  =  1.6  mm,  to  be  9.7°* 

Assume  that  the  shock  shape  is  stationary,  so  that  the  Mach  number  is  the 
same  along  the  shock  front.  At  each  point  on  it,  the  gas  is  deflected  toward 
the  shock,  as  shown  in  the  inset  of  Fig.  29.  The  overall  pattern  of  gas 
deflection  immediately  behind  the  shock  front  is  shown  by  arrows  which 
represent  the  direction  of  mass  flow.  Because  of  this  deflection,  compression 
occurs  at  C  behind  the  shock  front  with  positive  curvature,  and  expansion  occurs 
at  E  behind  the  shock  front  with  negative  curvature.  These  regions  C  are 
visible  in  Figs,  26  and  27=  This  compression  has  nothing  to  do  with  the  real- 
gas  rise  in  density  due  to  the  ionization  relaxation,  but  is  a  perfect-gas 
effect  at  the  curved  frozen  shock  front.  Compression  waves  form  and  converge 
to  weak  shocks  S.  Neighbouring  shocks  intersect,  forming  the  subsequent 
system  of  waves  that  appear  in  schlieren  pictures  as  thin  lines  (Fig.  26  and 
27)  and  in  interferograms  as  barely  discernible  fringe  shifts.  The  schlieren 
method  shows  up  these  transverse  waves  better  than  interferograms  (Fig.  22). 

Figure  26  also  shows  the  effect  that  the  ionization  relaxation  has  on  these 
transverse  waves.  After  the  shocks  intersect  twice  they  become  significantly 
steeper  and  remain  so  further  from  1^.  This  occurs  very  near  the  electron- 
cascade  front.  It  is  not  clear  what  is  happening  here,  since  the  nature  of  the 
transverse  waves  is  still  ambiguous. 

One  important  deduction  can  be  made,  however,  about  the  curved  shock 
front.  Regions  marked  C  in  Fig.  29  have  a  higher  than  average  pressure  due  to 
the  compression  caused  by  the  flow  deflection  introduced  by  the  curved  shock 
front.  Since  this  region  is  subsonic  relative  to  the  shock,  this  increased 
pressure  is  felt  at  the  shock  and  tends  to  increase  the  shock  speed;  thus 
lagging  regions  would  tend  to  catch  up.  Similarly,  leading  parts  tend  to  fall 
back.  This  seems  to  be  an  argument  for  stability  of  form,  but  the  influence  of 
the  relaxation  region  would  complicate  matters.  In  any  event,  the  part  that 
is  played  by  this  possible  mechanism  is  open,  since  stability  is,  in  fact, 
never  observed. 

Now  assume  that  the  shock  shape  evolves  as  a  standing  wave,  and  consider 
the  two  extremes  of  its  oscillation,  using  Fig.  26  as  an  example.  To  get  from 
1^  to  Ig,  the  shock  moved  to  the  left  with  an  average  Mach  number  Mg  (=l6.l). 

In  addition  to  this,  the  leading  and  lagging  parts  of  1^  travelled  different 
distances  in  the  same  time  At(=7.37  us)  to  reach  the  configuration  13.  This 
difference  in  distance  is  at  most  k  x  g  =  6.k  mm,  using  the  standing  wave  model 
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Eq.  3.4,  and  the  average  distance  between  1^  and  Ig  is  38.2  mm.  Thus,  one  point 
along  the  shock  front  could  attain  an  average  Mach  number  between  1^  and  Ig  as 
high  as  17*4,  while  another  point  could  have  one  as  low  as  l4.8.  Significant 
differences  in  frozen  and  equilibrium  pressure  and  temperature,  as  well  as  the 
Electron  number  density  and  the  relaxation  length,  occur  between  these  two 
Mach  numbers.  An  unsteady  two-dimensional  reacting-gas  shock  structure  would 
result.  This  rough  model  serves  only  to  gain  insight  into  the  complexity  of 
this  problem,  and  does  not  fully  exploit  the  standing-wave  model. 

Experimental  evidence  of  the  variation,  along  the  shock  front,  of  the 
relaxation  lenth  xE  can  be  observed  in  Fig.  21.  From  typical  profiles  of  ne 
through  the  shock  wave  (Ref.  2),  we  see  that  the  distance  between  the  shock 
front  and  the  place  where  dne/dx  is  a  maximum  is  nearly  equal  to  Xg.  The  region 
where  dne/dx  is  large  shows  up  on  this  schlieren  record  as  a  dark  band  D2.  This 
distance  between  the  shock  image  I,  and  the  band  is  least  when  1-^  lags  the 
mean  shock  position  and  greatest  when  1^  leads.  Since  we  know  that  Xg 
decreases  as  Mg  is  increased,  this  observation  implies  that  the  lagging  parts  of 
I,  have  been  travelling  faster  than  the  leading  parts.  This  implies,  in  turn, 
that  the  shock  image  had  a  greater  amplitude  at  a  previous  time.  Unfortunately, 
only  one  exposure  was  obtained  in  this  run  which  is  the  only  run  to  exhibit  a 
dark  band.  It  must  be  remembered,  however,  that  D2  is  the  result  of  the  geometry 
of  the  translational  shock  front  as  it  was  at  a  time  preceding  the  pulse  by 
the  induction  time  (the  time  taken  for  the  gas  to  reach  ionization  equilibrium 
after  passing  through  the  shock  front).  Thus,  such  an  extrapolation  could  be 
completely  unfounded.  However,  the  fact  that  D-^  follows  with  no  detectable 
shift  in  the  y-direction  supports  the  standing-wave  model.  Even  though  Fig.  21 
is  not  able  to  verify  the  standing-wave  model  regarding  physical  parameters 
such  as  xE,  it  can  help  stimulate  discussion  on  the  instability  phenomenon. 

Experimental  and  analytical  research  into  acoustic  stability  phenomena 
in  monatomic  plasmas  has  been  carried  out  almost  exclusively  on  electrically 
heated  plasmas,  where  high-energy  electrons  from  a  voltage  discharge  heat  the 
gas  to  ionization.  The  heating  is  done  by  the  net  transfer  of  energy  from  the 
electrons  to  the  atoms  by  elastic  and  inelastic  collisions,  where  Te  »  Ta. 

The  availability  of  hibh-voltage  sources  and  the  elegant  theoretical  treat¬ 
ments  possible  with  this  type  of  plasma  have  contributed  to  its  popularity. 

One  such  case  is  the  phenomenon  of  ion-acoustic  waves. 

On  the  other  hand,  the  plasma  behind  a  strong  shock  wave  in  argon  is  a 
thermal  plasma,  where  the  atom  temperature  is  suddenly  raised  significantly  in 
the  absence  of  any  free  electrons.  Ionization  occurs  initially  by  inelastic 
collisions  between  the  atoms  and  then  later  by  inelastic  collisions  between  the 
atoms  and  electrons.  Here  the  relation  Ta  >>  Te  holds  over  most  of  the  relax¬ 
ation  regions.  Appropriately,  the  only  literature  available  that  discusses 
thermal  plasmas  deals  with  the  more  specialized  topic  of  shock  structure  in 
strong  monatomic  gases  (Refs.  l,6,9»ll).  Except  for  papers  related  to  the 
presently-studied  phenomenon  (Refs.  2,10),  no  experiments  or  theoretical 
treatment  regarding  waves  in  a  thermally  heated  plasma  have  been  found. 

Morse  and  Ingard  (Ref.  7)  give  derivations  for  several  types  of  acoustic 
disturbances  in  the  electrically  heated  plasma  (Te  »  Ta),  further  dividing 
these  into  two  cases;  a  weakly  ionized  gas  (x  <<  l)  and  a  strongly  ionized  gas 
(X  ^  l).  The  former  implication  leads  to  a  dominately  neutral  gas  which 
receives  some  of  the  energy  that  the  electrons  draw  from  the  external  fields, 
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making  acoustic  feedback  and  amplification  possible.  The  latter  gives 
rise  to  "ion-acoustic"  waves  in  which  the  neutral  particles  may  be  ignored. 

In  each  case,  the  authors'  emphasis  is  on  the  effect  that  external  fields, 
such  as  a  strong  electric  field,  have  upon  the  charged  particles  in  the  plaBma. 
These  theoretical  considerations  and  the  present  situation  differ  in  three 
ways:  l)  we  have  here  Ta  »  Te  instead  of  Tg  »  T  ,  2)  we  have  a  intermediate 

value  for  the  degree  of  ionization  (0.03  <X  <  0.21),  so  that  neither  theory  is 
applicable,  and  3)  there  are  no  external  fields.  However,  since  the  net 
energy  exchange  between  species  of  different  temperatures  was  instrumental  in  the 
acoustic  disturbances  in  Ref.  7,  a  look  at  the  process  is  warranted  here. 

Expressions  for  the  net  energy  transfer  rate  between  electrons  and  atoms 
and  between  electrons  and  ions  in  elastic  and  inelastic  collisions  are  given 
in  Ref.  2.  The  energy  transfer  rate  for  elastic  collisions  between  electrons 
and  atoms  is  proportioned  to  ne,  na,  T  ,  and  Ta  -  Te.  The  corresponding  rate 
for  ions  instead  of  atoms  is  proportioned  to  ne  ,  Te  and  T&  -  Tg.  It  appears 
that  the  electron  number  density  ne  and  the  species  temperatures  Tg  and  T& 
probably  play  a  significant  role  in  determining  the  rate  of  energy  supplied  to 
the  plasma  distrurbance  which  creates  the  nonplanar  shock  front.  This  was 
realized  by  Bristow  (Ref.  l)  when  he  added  some  hydrogen  to  the  argon  to 
force  Ta  -  T  =  0,  and  hence  prevent  energy  transfer  to  the  hypothesized 
plasma  oscillation.  The  "resonant  cavity"  he  mentions  is  roughly  the  length 
of  the  relaxation  region  times  the  shock  tube  internal  crOss-sectional  area, 
and  this  cavity  follows  behind  the  shock  front  as  it  travels.  It  was  thought 
that  some  resonance  condition  existed,  depending  upon  the  volume  of  this 
"cavity"  and  the  energy  input  rate,  explaining  why  this  phenomenon  was  not 
observed  in  smaller  shock  tubes.  According  to  observations  in  this  study, 
discrete  resonance  conditions  seem  very  unlikely  owing  to  the  continuity  of 
njj  over  wide  ranges  of  shock  parameters.  In  other  words,  even  though  nh  has 
integer  values,  shock  with  a  particular  value  of  nh  occur  over  a  range  of  shock 
parameters.  Also,  evidence  of  a  translational  front  nonplanarity  with  an 
amplitude  comparable  (gk  from  Eq.  4,3)  to  those  in  Table  3  was  found  in  time- 
resolved  interferograms  in  Refs,  13  and  24.  The  value  of  nh  seems  to  be  1,  and 
the  amplitude  over  the  5  cm  channel  height  is  very  small,  about  1  mm.  The 
authors  make  no  mention  of  this,  most  likely  because  it  is  so  difficult  to 
notice  unless  looked  for.  However,  the  observation  of  the  effect  of  the  hydrogen 
impurity  still  stands,  suggesting  that  the  notion  of  some  plasma  oscillation, 
supported  by  energy  transfer  whose  rate  depends  upon  Ta  -  Tg,  is  valid. 

No  qualitative  description  can  be  given  of  the  plasma  oscillation  at 
this  time,  because  it  is  not  clear  what  type  of  oscillations:  l)  could  occur 
under  these  plasma  conditions  (values  of  ne,  x»  Ta,  Tg);  2)  would  be  consistent 
with  the  hydrodynamic  environment  of  the  shock  wave;  and  3)  give  rise  to  the 
nonplanar  translational  shock  front  as  observed  here  with  the  corresponding 
time  behaviour. 

Since  no  trend  in  the  amplitudes  g  of  the  instabilities  was  observed, 
attention  was  focussed  on  the  number  of  half-waves  njj  in  an  attempt  to  link  a 
characteristic  of  the  instability  to  the  physical  parameters  of  the  shock, 
most  notably  ne  and  T  .  As  mentioned  in  Sec.  2.1.2,  because  the  initial 
temperature  T,  is  fixed  (room  temperature),  the  shock  wave  produced  in  argon  is 
completely  determined  by  the  two  independent  pressures  p1  and  p^  (from  Pj,hx). 

From  these  two  values  one  can  calculate  the  shock  Mach  number  Mg  (Eq.  2.3),  and 
from  the  conservation  equations  and  Saha  equation  all  the  flow  and  thermodynamic 


variables  n  ,  n&,  x»  Tg,  etc°  of  this  complete  set 
choose  any  two  as  the  independent  ones.  Thus,  to  determine 
determine  p-y  and  p^,  exactly  two  of  these  variables  must  be 
practise  this  is  how  a  run  was  planned. 


of  variables  we  may 
a  run,  that  is, 
chosen.  In 


In  view  of  the  importance  of  ne  and  T  -  Te  in  the  expressions  for  the 
energy  transfer  associated  with  possible  plasma  oscillations,  each  run  was 
planned  by  choosing  the  maximum  electron  number  density  ne,E  at  quasi-equilibrium 
and  the  frozen  temperature  T2(fr)  the  translational  shock  front.  Although 

the  value  of  ne  is  zero  at  the  shock  front  and  rises  to  ne,E,  it  was  felt  that 
the  magnitude  of  ng,E  would  give  an  indication  of  ng  throughout  the  region  where 
Ta  >  Tg.  For  simplicity,  Ta  is  considered  instead  of  T&  -  Tgi  To(fr)  is  the 
atom  temperature  immediately  behind  the  shock  front  for  most  of  the  relaxation 
length  (Ref.  2),  so  T&  may  be  approximated  by  Tg ( f r ) •  Since  Tp(fr)  Just  the 
perfect-gas  temperature  behind  the  shock  front,  it  depends  only  on  Ms  since  T^ 
is  practically  fixed.  This  justifies  choosing  ng,E  and  Ms  as  the  independent 
parameters  of  the  shock.  The  value  of  pj  was  found  from  thermodynamic  tables, 
and  pjj  was  obtained  from  Eq-  2.3-  The  experiment  was  then  performed,  and  the 
observed  M  and  initial  pressure  p^  were  measured  (with  accuracies  given  in 
Sec.  2.1.1+}.  To  add  the  experimental  data  to  Fig.  30,  a  thermodynamic  routine 
was  used  to  calculate  ne,g  from  these,  and  a  symbol  representing  the  value  of 
njj  observed  from  the  schiieren  picture  was  placed  at  the  point  given  by  Mg  and 
ne,E  (all  three  for  each  run  are  given  in  Tables  1  and  2),  In  this  way  trends 
regarding  the  values  of  ng,E  and  Ms  for  one  value  of  would  be  seen.  Figure  lk 
serves  to  give  the  region  containing  values  of  ng,E  and  Ms  attainable, 
consistent  with  the  constraints  of  the  driver  ip^  £  1+43  atm  (6500  psi)],  the 
shock  detectors  (p^  >  1  torr),  and  the  test-section  windows  lp2  <  6.8  atm  (100 
pse)}.  As  well,  experiments  were  designed  to  get  sample  experiments  from  as 
many  parts  of  the  ng,E  -  Ms  plane  as  possible,  ir.  order  to  obtain  comprehensive 
data  on  the  phenomenon.  It  must  be  emphasized  here  that  the  position  of  a 
plotted  point  in  Fig.  30  has  significance  cnly  with  regard  to  the  value  of  nh. 
Thus,  curves  formed  by  points  with  different  values  of  ny,  have  no  meaning, 
whereas  a  distribution  of  points  with  the  same  n^  suggests  a  relationship 
amongst  ne,E,  Ms  and  nE , 


Some  trends  can  immediately  be  observed  from  Fig,  30.  For  some  constant 
Ms,  nyj  increases  with  ne,E  almost  in  proportion.  Second,  for  some  constant 
ng,E,  n^  increases  with  Mg.  Remember  that  the  uncertainty  of  n^  is  ±1. 

Despite  the  scatter  of  points  with  the  same  value  of  n^,  a  general 
relationship  amongst  nyj,  ne,E,  and  can  also  be  seen.  Its  most  conspicuous 
manifestation  is  the  accumulation  oi'ny,  =  6  in  the  upper-right  corner  of  the 
concentration  of  nh  -  2  in  the  lower-left.  These  distributions  are  reflected 
in  the  diagonal  lines,  one  for  each  nh,  drawn  such  that  an  equal  number  of 
points  with  that  value  of  n^  fall  on  either  side.  A  regularly  spaced  pro¬ 
gression  from  nh  ~  1  m  the  lcwer-left  to  nh  -  6  at  the  upper-right  is 
indicated,  except  for  n^  =  it.  Other  than  this  anomaly,  it  seems  that  the 
value  of  n^  can  be  predicted  for  a  given  experiment  within  ±1  in  this  facility 
given  the  value  of  ng,E  and  Mc  If  the  p^  -  grid  from  Fig,  lit  is  superposed 
on  Fig  30,  nh  can  be  predicted  within  il  in  term?  of  the  initial  experimental 
pressures  p-y  and  p(^, 


Formulating  this  prediction  in  terms  of  the  straight  lines  in  Fig.  30 
(analogous  to  nh  "contours")  would  result  in  equations  with  no  physical 
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foundation.  This  would  be  remedied  by  determining:  l)  the  physical  process 
which  gives  rise  to  the  instabilities:  2)  the  dependence  of  this  process  upon 
the  various  parameters  such  as  ne,  n  ,  x»  T  ,  and  Ta;  and  3)  the  link 

between  the  instability  characteristic  n^  and  the  relevant  parameters.  No 

solution  to  these  steps  is  offered  here.  To  investigate  n^  as  a  function  of 
these  variables  several  combinations  of  ne,  na,  Te,  and  Ta  were  considered 
keeping  in  mind  their  roles  in  the  previously-mentioned  energy  transfer  rates. 
Comparing  the  values  of  these  expressions  over  runs  with  the  same  observed  n^ 
met  with  less  success  than  the  simple  graphical  predictor  of  Fig.  30.  However, 
the  possibilities  are  by  no  means  exhausted. 

As  a  predictor  of  the  geometry  of  the  nonplanarity  of  the  shock  front  from 
given  ne,g  and  M  ,  Fig.  30  applies  only  to  this  shock  tube  facility.  A  more 
universal  indicator  would  be  the  wavelength  of  the  sinusoidally  curved  shock- 
front  disturbance,  given  by  X  =  2L/njj,  where  L  is  the  channel  height.  The 

physical  meaning  of  X  is  more  amenable  to  nonintegral  values  of  n^,  which 

woule  contradict  the  wall  condition  but  would  be  no  more  troublesome  than  the 
travelling-wave  model.  If  X  is  used  instead  of  n^  then  a  quantitative  com¬ 
parison  could  be  made  with  other  shock  tubes.  Thus  it  seems  reasonable  that 
in  a  smaller  shock  tube  say  with,  a  channel  height  of  L  =  6  cm  one  would 
observe  n^  =  2  for  a  shock  wave  that  in  this  facility  would  appear  as  nfa  =  6, 
the  wavelength  X  of  the  disturbance  being  the  same.  The  values  of  X  are 
listed  for  each  run  in  Table  3. 

In  Table  3  are  listed  the  apparent  periods  of  oscillation  of  the  shock- 
front  instability.  Each  travelling  wave  has  a  period  of  about  15  ys.  Three 
standing  waves  have  periods  of  the  magnitude  50  -  90  ys  and  lower  n.,  whereas 
the  other  two  have  periods  also  near  15  ys  with  n^  =  6.  Note  that  the  longest 
periods  are  for  runs  with  the  lowest  ne,E,  whereas  the  shortest  period  occurs 
for  the  run  with  the  highest  ne,E  and  a  measureable  period.  This  observation 
is  consistent  with  the  notion  that  a  higher  electron  number  density  in  the 
reacting  region  of  the  shock  would  increase  the  rate  of  energy  transfer  to 
the  electrons.  No  such  correlation  can  be  seen  using  the  amplitude. 

Conclusions: 

1)  It  is  possible  to  obtain  multiple  short-duration  monochromatic  light 
pulses  from  a  ruby  laser  in  a  combined  giant-pulse/normal-lasing  mode 
with  a  flash-lamp  pumping  time  of  about  500  ys,  the  first  two  pulses 
consistently  being  about  7  ys  apart. 

2)  Using  a  very  sensitive  schlieren  system  with  a  dark  background,  this 
light  source  can  be  used  to  photograph  a  high-speed  low-density 
ionizing  shock  wave  in  argon  at  two  discrete  times.  This  apparatus 
was  used  to  investigate  the  time -dependent  behaviour  of  the  unstable 
translational  shock  front  observed  in  ionizing  shocks  in  argon. 

3)  The  shock-front  nonplanarities ,  which  consistently  appear  as  smooth 
sinusoidal  curves,  were  observed  to  change  in  time  either  like  a 
travelling  wave  or  like  a  standing  wave,  with  periods  of  oscillation 
decreasing  and  number  of  wavelengths  increasing  with  increasing  ne,E. 


U)  A  comprehensive  series  of  experiments  led  to  a  graphic  predictor 

designed  to  give  within  ±1  the  number  of  half-waves  in  the  nonplanar 
front  of  a  given  shock  wave  in  argon.  In  terms  of  wavelength  of  the 
sinusoidal  disturbances,  the  predictor  might  be  applied  to  other 
facilities. 

5)  An  extensive  system  of  transverse  waves  was  observed  behind  some 

of  the  curved  shock  fronts.  A  possible  hydrodynamic  explanation  is 
given  which  relates  these  waves  directly  to  the  geometry  of  the 
shock  front. 

6)  Additional  analyses  are  required  in  order  to  explain  this  phenomenon 
in  terms  of  the  ionizing  shock  properties  and  the  experimental  data 
obtained  in  the  present  investigation. 
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TABLE  3 


GEOMETRICAL  PARAMETERS  OF  THE 
TIME-DEPENDENT  SHOCK-WAVE  NONPLANARITIES 


Run  # 
from 
Table  2 

t 

P 

:  (ys) 

n 

s 

At 

(ys) 

max  g  n^ 

(mm) 

\ 

(cm) 

T 

(ys) 

Wave 

Type 

1 

1*82 

2 

7.58 

0.7 

2 

17.8 

2 

1+79 

1* 

6.86 

1.8 

1 

35.6 

1+8. 

S 

3 

5U8 

1 

1.5 

2 

11.8 

1* 

1*92 

3 

6.78 

1.5 

2 

17.8 

S  (or  T) 

5 

1*90 

1 

0.8 

3 

11.7 

6 

536 

1 

1.5 

1 

35.6 

7 

1*90 

3 

6.57 

3.2 

2 

17.8 

59. 

S 

8 

50U 

1 

1.6 

3 

11.8 

9 

1+1+3 

7 

6.99 

1-.0 

2 

17.8 

10 

511 

2 

12. 

1.0 

2 

17.8 

S 

11 

1+90 

3 

7.57 

1.7 

2 

17.8 

12 

510 

2 

7.80 

1.6 

5 

17.1 

13 

1*1*6 

3 

6.61* 

3.3 

2 

17.8 

S 

lit 

1+93 

2 

7.87 

1.5 

1* 

8.9 

S? 

15 

1+77 

3 

6.9I+ 

2.2 

3 

11.8 

88.? 

S 

16 

1+77 

2 

7.37 

1.2 

3 

11.8 

16.3 

T 

17 

1*62 

2 

7.1+2 

1.1* 

5 

7.1 

S? 

18 

501 

2 

8.1 

1.0 

6 

5.9 

16. 

S 

19 

1*83 

1* 

7.7l* 

1.1+ 

1* 

8.9 

20 

1*96 

3 

7.37 

1.6 

6 

5.9 

15. 

S 

21 

1*67 

5 

7.1+1+ 

1.8 

5 

7.1 

15.3 

T 

22 

512 

1* 

7.98 

1.7 

6 

5.9 

1I+.6 

T 

23 

515 

1 

1.2 

6 

5.9 

*P- 

laser  flash-lamp  pumping  time 

»  nS  “ 

number 

of  visible  shock  images; 

At  - 

time  between 

P]_  and 

p2;  nh  - 

number 

of  half-wavelengths 

» 

T  -  possible  period  of  sinusoidal  non-planarity; 

Wave  Type  -  standing  (S)  or  travelling  (T)  sinusoidal  wave. 


FIG.  1.  COORDINATE  SYSTEMS  IN  THE  SHOCK  TUBE.  x,y,z:  LABORATORY 

COORDINATES.  x'.y'.z':  SHOCK-FIXED  COORDINATES. 


INTERFEROGRAM  OF  THE  ONE  AND  TWO-DIMENSIONAL  SHOCK  FRONTS  IN 
A  SHOCK  TUBE  (REF.  18).  S  -  STRAIGHT  ONE-DIMENSIONAL  SHOCK 
WAVE  (FRINGE  SHIFT  IS  NEARLY  AN  INTEGER).  C  -  CURVED  TWO- 
DIMENSIONAL  SHOCK  WAVE. 


FIG.  2. 


3  IilTERFER  OGHAM  OF  A  THREE-DIMENSIONAL  SHOCK  FRONT  MOVING  LEFT 
IN  IONIZING  KRYPTON.  NOTE  THE  'TWISTED  RIBBON'  APPEARANCE 
OF  THE  TRANSLATIONAL  SHOCK  FRONT.  THE  SHOCK  SURFACE  DEPENDS 
ON  y',  AND  z'  (REF.  lc)  . 


FIG.  4  CONGRUENT  FRINGES  PRODUCED  BY  PLANAR  SHOCK  FRONT  IN  IONIZING  ARGON 
MOVING  LEFT  (REF.  10). 


riON-PIANAK  SHOCK  :  RONT  IN  IONIZING  ARGON  MOVING  LEFT 
SHOWING  SINUSOIDAL  OSCILLATIONS  (RUN  B5  IN  TABLE  l)  . 
Ma  -  lb .01,  Vi  =  l.-W  TOKR,  Tx  =  295. 3’K. 


6  PLANAR  SHOCK  WAVE  AND  INDUCED  PLANAR  FLOW  IN  DISSOCIATING 
NITROGEN  (REF.  l)  .  Ms  =  16.6,  px  =  2.26  TORR,  Tj,  =  295-6°K. 


FIG.  7(a)  ATOM  TEMPERATURE  Ta  AND  ELECTRON  TEMPERATURE  Tg  BEHIND  A 
STRONG  INCIDENT  SHOCK  WAVE  IN  ARGON  FOR  Ms  =  18.0,  pp  = 

3  TORR,  Tl  =  296. 5°K  (REF.  23). 


FIG.  7(b)  RATIOS  Ta/Ti  AND  Te/Ti,  AND  DEGREE  OF  IONIZATION  BEHIND  A 
STRONG  INCIDENT  SHOCK  WAVE  IN  HYDROGEN  FOR  Ms  =  25,  pi  = 

3  TORR  (REF.  25). 


VS 


x ; . 


il 


t  0! R-.  IHrKRFEWV'.KAM'.  .'.t  APART  FROM  IMA,;K  CONVERTER  CAMERA  O'  THE 
SAME  .'H'lCK  WAVE  MOVIN' :  LEFT  IFTTO  AROOIJ.  Tiff:  TOT  l:-AMK  WAV  EXPOSED 
i'lKrST,  TiK  BOTTOM  FRAME  LAST.  EACH  CTUCTLAR  NE  iATIVF:  IMA  IE  WAT 
1.  MM  IN  HIA. 


(a)  Ma  -  1(. .L’,  px  =  5.15  TORR,  At  = 

(b)  Ft,  =  10.35,  P1  =  5.U  TORR,  Ft 

in  each  Kim  tiif:  shock  front  can  he  seen 

TOP  TO  THE  BOTTOM  FRAME. 


US  . 

7  us. 

TO  M HVE  I£?T  KHO M  THE 


Measured  Initial  Mixture  Pressure  (H2  ,02 ,  He  )  =  423  psi  (Heise  Gauge) 
Calculated  Diaphragm  Bursting  Pressure  =  2452  psi 

*|  |—  5  msec 


Final  Combustion  Pressure 


Initial  Mixture  Pressure 

Original  Base  Line  — 
P»0 


1 

Arrival  of  Rarefaction 
Wave  after  Diaphragm 
Burst 

1 

500  psi 


a)  DRIVER  PRESSURE  at  CLOSED  END 

Measured  Shock  Mach  Number  =  15.0  at  Station  M 
Measured  Initial  Pressure  Po  =  6.l2  torr  (0.12  psi) 


Predicted  Equilibrium 
Pressure  *  38.0  psi 
(1965  torr) 

Predicted  Frozen 
Pressure  =  33.3  psi 
( 1722  torr  ) 


J 


10  psi 


10  yu.  sec 


b)  SHOCK- WAVE  PRESSURE  at  STATION  M 


FIG.  10.  PIEZO  PRESSURE  HISTORIES  IN  COMBUSTION  DRIVER  a)  AND  IN  TEST 
SECTION  (b). 


FIG.  11(b)  DISTANCE-TIME  CURVE  OF  THE  INCIDENT  SHOCK-WAVE  AS  IT  MOVES  ALON3 
THE  SHOCK  TUBE.  THE  SIX  POINTS  CORRESPOND  TO  STATIONS  D,  F,  G, 

H,  I,  AND  J  OF  FIG.  11(a).  TIMES  ARE  HEASURED  FROM  D.  THE  SOLID 
LINE  REPRESENTS  THE  BEST  FIT  ACROSS  THESE  MEASUREMENTS.  THE 
VELOCITY  OF  THE  INCIDENT  SHOCK-WAVE  AT  ANY  POINT  IS  dx/dt. 


FIG.  14.  VALUES  OF  np,r  vs  Hs  ATTAINABLE  IN  THE  UTIAS  SHOCK  TUBE  FOR  ARGON 


OUTPUT  OF  PHOTODIODE  EXPOSED  TO  LASER  LIGHT  WHEN 
VERTICAL  -  RELATIVE  INTENSITY,  HORIZONTAL  -  TIME 
Px  TO  Pj  ARE  OUTPUT  SIGNALS. 
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Fid.  20(a)  MULTIPLE- EXPOSURE  3CHLIEFEN  PICTURE  OF  AM  IONIZING  SHOCK  WAVE  MO'/; 

LEFT  INTO  ARGOII.  Ms  *  11.5,  Pi  =  9.98  TORR,  Tl  =  299. 2°K,  tp  ” 

US  (RUN  NO.  2).  TWO  USABLE  SHOCK  IMAGES  APPEAR  AT  I !  AND  I2;  TWO 
FAIIfT  ONES  AT  I3  AND  11,  .  ONLY  THE  TRANSLATIONAL  FRONT  IS  IMAGED 
HERE,  APPEARING  AS  A  CURVED  'VERTICAL  LINE. 


FIG.  20(b)  SCHEMATIC  DRAWING  OF  THE  CURVED  SHOCK  IMAGE  I,  AT  LEFT,  SHOWING 
THE  LINES  £  AND  i-[  WHICH  ARE  PARALLEL  TO  (THE  Y-DIRECTION) , 
THE  SHOCK  POSITION  MIDWAY  BETWEEN  THEM,  AND  THE  AMPLITUDE  g 
OF  THE  DISTURBED  SHOCK  FRONT,  EXAGGERATED  HERE  FOR  THE  PURPOSE 
OF  ILLUSTRATION.  IN  (a)  2g  =  3.7  mm.  Ii  CONSISTS  OF  THE  BASIC 
GEOMETRIC  UNIT,  THE  COSINE  HALF-WAVE,  WHICH  IS  PERPENDICULAR  TO 
THE  TOP  AND  BOTTOM  WALLS.  NOTE  THAT  Ix  CROSSES  li  JUST  ONCE. 


W1  -  WAVE  SYSTEM  BEHIND  THE  CURVED  SHOCK  FRONT,  D*  -  DARK 
BAND  DUE  TO  LARGE  dne/dx  AT  ELECTRON  CASCADE  FRONT.  NOTE 
THAT  THIS  BAND  VARIES  REGULARLY  IN  DISTANCE  FROM  It  ACCORDING 
TO  THE  LATTER'S  SHAPE. 


MULTIPLE-EXPOSURE  3CHLIEREM  PETURE  OF  AN  IONIZING  SHOCK  WAVE 
MOVING  LEFT  III  ARGON.  Ms  =  1** . 5 ,  Pi  =  3.18  TORR,  Tj  =  299-0°K 
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FIG.  24(b)  SCHEMATIC  DRAWING  OF  THE  IMPORTANT  FEATURES  OF  24(a). 

Il  TO  I3  -  SCHLIEREN  IMAGE  OF  THE  TRANSLATIONAL  SHOCK 
FRONT,  ij,  TO  ^3  -  MEAN  SHOCK  POSITIONS,  wy  -  VERTICAL 
REFERENCE  WIRES.  ALL  THREE  SHOCK  IMAGES,  ESPECIALLY 
II,  EXHIBIT  THREE-DIMENSIONAL  EFFECTS. 


FIG.  25(b)  SCHEMATIC  DRAWING  OF  THE  IMPORTANT  FEATURES  OF  25(a). 

Il  -  SCHLIEREN  IMAGE  OF  THE  TRANSLATIONAL  SHOCK  FRONT. 

Wt_  -  SYSTEM  OF  TRANSVERSE  WAVES  BEHIND  I]_ .  Di  -  ELECTRON 
CASCADE  FRONT. 
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FIG.  29(a)  FLOW  DEFLECTIONS  BEHIND  SINUSOIDALLY  CURVED  TRANSLATIONAL  IONIZING 

SHOCK  FRONT  AND  SUBSEQUENT  TRANSVERSE  WAVE  SYSTEM.  GAS  ENTERS  SHOCK 
FRONT  j'ROM  LEFT  MASS  FLOW  REPRESENTED  BY  ARROWS  IS  DEFLECTED 
ACCORDING  TO  ANGLE  OF  INCIDENCE  p;  DIVERGING  ARROWS  RESULT  IN 
EXPANSION  REGIONS  E;  CONVERGING  ARROWS  GIVE  RISE  TO  COMPRESSION 
REGIONS  C,  WHICH  CONVERGE  TO  FORM  TRANSVERSE  SHOCKS  S;  SHOCKS  NEAR 
THE  WALL  REFLECT,  MAINTAINING  PATTERN.  IF  SHOCKS  ARE  OF  EQUAL 
STRENGTH  NO  CONTACT  SURFACE  RESULTS;  OTHERWISE  CONTRACT  SURFACES 
WILL  FORM. 

(b)  FLOW  ENTERING  AN  OBLIQUE  SHOCK  IS  DEFLECTED  TOWARD  THE  SHOCK: 

NORMAL  VELOCITY  ACROSS  SHOCK  WAVE  IS  REDUCED;  PARALLEL  COMPONENT 
REMAINS  UNCHANGED. 
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